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ACTIVE FEED ARRAY COMPENSATION FOR REFLECTOR ANTENNA 
SURFACE DISTORTIONS 

Roberto 0. Acosta 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 

SUMMARY 

The feasibility of electromagnetic compensation for reflector 
antenna surface distortions (non-random) Is Investigated. The 
performance characteristics (gain, sldelobe levels, etc.) of large 
satellite communication reflector antenna systems degrade as the 
reflector surface distorts, mainly due to thermal effects from a solar 
radiation. The technique developed In this dissertation can be used 
to maintain the antenna bores Ight directivity and sldelobe level 
Independent of thermal effects on the reflector surface. With the 
advent of monolithic microwave Integrated circuits (MMIC), a greater 
flexibility In array-fed reflector antenna systems can be achieved. 
MMIC arrays provide Independent control of amplitude and phase for 
each of the many radiating elements In the feed array. By assuming a 
known surface distortion profile, a simulation study Is carried out to 
examine the antenna performance as a function of feed array size and 
number of elements. Results Indicate that the compensation technique 
can effectively control boreslght directivity (within 1-3 dB of the 
undistorted) and sldelobe level (within 1-5 dB of the undistorted) 
under peak surface distortion In the order of tenth of a wavelength. 
Furthermore, the computational ease of the compensation technique 
allows It to be Implemented adaptively In a large satellite reflector 
antenna system. 
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CHAPTER 1 
INTRODUCTION 

Future space communication multibeam antenna system will utilize 
large reflector surfaces. Maintaining the surface accurracy of this 
large reflector antennas is important, particularly when high gain 
(higher than 50 dB) and low sidelobes (lower than -30 dB) are 
required. Reflector surface errors are generally classified into the 
following types: 

(i) Random surface error in the order of several mils due to 
manufacturing imperfection (Refs. 1-3). 

(ii) Large-scaled surface error in the order of inches due to 
thermal, gravitational and other effects (Ref. 4). 

It is the large-scale surface error that is of interest in this 
dissertation work. Methods for reducing this type of surface error 
include the use of better antenna supporting structure, better 
thermally-insensitive material, and mechanically tunable surfaces. A 
different approach is to compensate for the degraded antenna performance 
due to surface distortions by using an active phased array feed. The 
later approach is becoming increasingly attractive due to recent 
advances In monolithic microwave integrated circuits (MMIC). MMIC 
arrays provide independent control of amplitude and phase for each of 
the many radiating elements in the feed array. Surface error 
reduction methods (mechanical/material) and surface error 


compensation techniques (electromagneti cally) are complementary, each 
approach can further improve the antenna performance independently of 
the other. 

This dissertation investigates surface error compensation 
techniques for large scale surface distortions. In chapter 2, 
techniques for calculating the reflector antenna secondary pattern are 
presented. A special emphasis is placed on the aperture integration 
and the physical optics methods. These techniques are necessary tools 
in evaluating the reflector antenna performance (directivity, sidelobe 
levels, etc.). . 

In chapter 3, a method for uniquely representing the surface 
errors of a distorted reflector antenna is developed. The distorted 
reflector surface (numerically specified) is separated into two 
components: an undistorted surface component (represented by a best 

fit paraboloid) and a surface error component (expanded into a fourier 
series). Correlation between surface error spectral components and 
antenna radiation performance is examined. The result is compared 
with the ones obtained by using a spline interpolation algorithm. 

Chapter 4 presents one of the two approaches for calculating the 
compensating feed array complex excitations, namely the indirect 
conjugate field matching (ICFM) approach. It is "indirect" because 
the compensating feed array excitations are determined by using the 
transmitting mode and the reciprocity theorem. The design of the 
compensating feed array (array size, number of elements, etc.) and its 
relation to the surface distortion are examined. 
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Chapter 5 presents the second compensation approach based on the 
receiving mode, namely the direct conjugate field matching (DCFM) 
approach. It is "direct" because the compensating array excitations 
are determined by integrating directly the induced surface currents on 
the reflector under a receiving condition. In contrast to the first 
approach, the DCFM allows the reflector be illuminated by an incident 
plane wave with a taped amplitude. The level of taper can effectively 
control the sidelobe level of the compensated antenna pattern. A case 
study is presented for comparing results between the DCFM and ICFM 
techniques. 

Finally, concluding remarks are given and some future work 
evolving from this dissertation research are discussed in chapter 6. 
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CHAPTER 2 


REFLECTOR ANTENNA SECONDARY PATTERN COMPUTATION 
2.1 Introduction 

Reflector antennas are widely used in communication satellite 
systems because they provide high gain at low cost. In analysing 
reflector antennas the computation of the secondary pattern is the 
main concern. The basic techniques for calculating the scattered 
vector fields from a reflector surface has been extensively developed 
in the open literature (Refs. 5 to 25). The best known of these 
techniques are the ones based upon aperture integration and physical 
optics methods. 

(i) Aperture Integration Method (AI, Fig. 1(a)). The fields 
on the aperture plane S ap are first calculated by tracing a 
reflected ray P] 0 r P 2 using geometrical optics (GO) theory and an 
edge diffracted ray Pi 0^ P 2 using geometrical theory of diffraction 
(GTD). The fields over the aperture plane 5Z ap are then integrated 
via a fast fourier transform algorithm (FFT) to obtain the scattered 
vector fields. 

(ii) Physical Optics Method (PO, Fig.. 1(b)). The induced 

currents on the reflector surface is approximated by 2 • r? x ifi . 

Where H 1 ! is the incident magnetic field intensity on the surface, 
and ^ is the surface normal. An integration of this current over 
the curved reflector surface gives the scattered vector fields. 






The exact solution of the scattered vector field can be expaned 
in a high-frequency asymptotic series (Refs. 26 to 27), namely, 
f?S(r) ~ ^-JKS(r) ( J o(r) / <-l /2 ^ (r) + H ^ 2 (r > + - • • • > for < * » 

( 2 . 1 ) 

Where K is the free space wavenumber. 

The two pattern computational methods described above recover terms in 
(2.1) to various degrees, as summarized in Table 1. 


TABLE 1. - ACCURACY AND LIMITATIONS OF PO AND AI 


Methods 

Accuracy 

Limitation 

AI (with GO rays 
and no edge 
diffracted rays) 

Recover X 0 

Not accurate for 
wide-angle sidelobes 

PO 

Recover 7£ 0 and 
partial A] 

AI (with GO rays 
and edge dif- 
fracted rays) 

Recover A 0 and Tt] 

Infinite field at 
caustics 


A detailed description of the computation of the incident field 
and the power radiated by the feed source is presented in section 2.3. 
The AI and PO methods are described in sections 2.4 and 2.7 
respectively. A numerical discussion is presented in section 2.8. 

2.2 Description of Problem 

The goal of reflector antenna secondary pattern computation can 
be described as follows: Given the feed element characteristics 

(location, polarization and element pattern) and reflector surface 
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characteristics (x,y,z, its first and second derivatives), it is 
required to compute the scattered vector field at an arbitrary 
observation point in the far-field zone. The geometry of the problem 
under consideration is depicted in Fig. 2. 

The reflector is shown in the main coordinate system (x,y,z). 
The feed coordinate system (xf,yf, Zf) facilitates the description of 
the feed polarization and radiation characteristics. This coordinate 
system is related to the main coordinate system by a set of Eulerian 
angles (Ref. 29). The time factor is suppressed throughout. 

It is assumed that the feed source has a well defined "phased center" 
at P] with coordinates ( x -j ,y i , z-j ) , and radiates a spherical wave 
denoted by (T^,.^). If the feed source is an array, each element in 
the array is considered separately. The secondary pattern of the 
array-fed reflector antenna is then obtained by vectorial ly adding the 
individual element scattered fields. The conducting reflector surface 
may be described by an analytical equation or by a set of points. For 
the numerically specified case, an interpolating technique (Refs. 30 to 
32) should be used for obtaining the required surface derivatives 
(first and second derivatives) . 

2 . 3 Incident Field Computation 

The radiated electric and magnetic fields from the feed source 
have the following asymptotic form 

-Jl ( } „ e_'^ Kr f ? (0 f , <|> f ) (2.2a) 

f r f 


^(r f ) ~ 


x tUr*) 


Z 


(2.2b) 



SCATTERED ELECTROMAGNETIC 



2 Geometry for calculating the scattered electromagnetic vector 
fields from a reflector antenna 
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Where (rf, ©f, <j>f) are spherical coordinates of an arbitrary 
observation point in the feed coordinate system, Z is the wave 
impedance, (6f, <j>f) is a vector function describing the feed source 
radiation pattern, K = 2.ir/\ is the wavenumber, X is the wavelength, 
and rf is the distance from the source to an arbitrary point in the 
reflector surface. The vector function (6f, <|>f) in Eq. (2.2a) can 
be approximated (Ref. 23) by Eq. (2.3). 

^(6f. <j>^) = ©^ UE(6^) (ae^ cos <j>^ + b sin 4y) + 

9f UH(0 f ) (b cos <jy - ae^ sin <j>^) (2.3) 

where UE (6f) is the feed E-plane pattern and UH (Of) is the 
corresponding H-plane pattern. a,b, and <|» are the feed 
polarization parameters: 


TABLE 2. - FEED POLARIZATION PARAMETERS 


Polarization type 

a 

b 

t 

Linear X 

i 

0 

O 

Linear Y 

0 

1 

0 

Right-hand circular (RHCP) 

0.707 

0.707 

0 

O 

o> 

Left-hand circular (LHCP) 

0.707 

0.707 

o 

O 

cn 


Typically the element pattern can approximated by a (cos (©f ) ) c l 
function that is. 
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UE (6 f ) = (cos (6f)>q E (2.4a) 

UH (6 f ) = (cos (Sf))^ (2.4b) 

If Eqs. (2.4a) and (b) are used to represent the feed 
radiation pattern, the power radiated by a single feed source is given 
as 


p (qE + qH + 1) 

rad ~ 60(2qE + 1)(2qH +1) 

For a feed array with M elements (Ref. 33) the power radiated is 
given by Eq. (2.6) . 


(2.5) 


P 

' rad 


n m 

■ I I 

m=l n=l 


*m*n \n 


( 2 . 6 ) 


where I m is the m-th element complex excitation, A is a square 
matrix. The derivation for the matrix A is included in appendix A. 
The radiated power of an antenna is an important quantity, since it is 
used in the computation of its directivity (section 2.8). 

Figure 3 illustrates the Eulerian angles (G1.G2.G3) which 
determine the rotational relationship between the main coordinate 
system (x,y,z) and the feed coordinate system (xf.yf.Zf). The 
definition of these angles are as follows. The angle G1 describes a 
counterclockwise (ccw) rotation about the z axis which brings the x 
axis to the x" axis aligned with the line of nodes (line of 
intersection between x-y and xf-yf planes), angle G2 defines a 
rotation about the line of nodes in a ccw sense as indicated so that 
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this brings the z axis to Zf, and angle G3 is another rotation about 
the Zf axis and aligns the x" axis with the xf axis in a ccw 
sense. 


In general the feed coordinate system orientation is expressed 
relative to the main coordinate system by a set of three orthogonal 

, v-\ ^ ^ 


unit vectors (xf,yf,zf). Let these unit vectors be expressed by: 

Xf = axl x + ax2 y + ax3 z 

yf = ayl x* + ay2 y“ + ay3 z~ 

Zf = azl T + az2 y* + az3 IT 

Then, the Eulerian angles are given by: 

azl 


G1 = tan' 

G2 = tan' 
1 / 

G3 = tan 


-az2 




ax3 


It can be shown that 
where 


B = 


yazl 

x ax2 

- ax 

z?)T = 

B(x y 

‘z)T 

bll 

bl 2 

bl 3 

b21 

b22 

b23 

b31 

b32 

b33 


(2.7a) 

(2,7b) 

(2.7c) 

(2.8a) 

(2.8b) 

(2.8c) 

(2.9) 


( 2 . 10 ) 


in which the elements of the Eulerian matrix B are defined as 
fol lows : 


bn = cos G1 cos G3 - sin G1 cos G2 sin G3 (2.11a) 

b 1 2 = s ^ n G1 cos G3 + cos G1 cos G2 sin G3 


(2.11b) 
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b 13 = 

sin 

G2 

sin G3 



(2. 

11c) 

b 21 = 

- cos G1 sin 

G3 

- s 

in G1 cos 

G2 

cos G3 

(2. 

lid) 

b 2 2 = 

- cos G1 cos 

G2 

cos 

G3 - sin 

G1 

sin G3 

(2. 

lie) 


b 2 3 = 

sin 

G2 

cos G3 



(2. 

,llf) 


b 31 = 

sin 

G1 

sin G3 



(2. 

1 1 g) 


b 32 = - 

cos G1 

sin G2 



(2. 

. 1 1 h) 


b 33 = 

cos 

G2 



(2. 

ni) 


x fs 


"x s - X]‘ 

Vfs 

— 

= B 

y s - yi 

-Zfs- 

’ 

• z s - z l- 


Thus, a point on the reflector surface with coordinates P s : 

(x s ,y s ,z s ) in the main coordinate system is transformed into the feed 
coordinate system through the transformation. 


( 2 . 12 ) 


where (xi.yj.Z]) is the origin of feed coordinate system indicating 
its transformation with respect to the main coordinate system. The 
corresponding spherical coordinates for this point is given by 

Zf s'-) (2.13a) 

(2.13b) 
(2.13c) 

The incident electric and magnetic fields are obtained from 
Eqs. (2.2a), (b) and (2.13a) to (c). Notice that (rf) and iturf) 
are in the feed coordinate system. Let the spherical vector components 
of l:Urf) and "ph(rf) be denoted by 


\/(Xfs 2 + 

yfs 2 

-1 

( z f S 

y = cos 

\ r f 

^ -1 

n f S 

y = tan 

fe 



1 4 


^(rf) = (E r f, Eef , E^f ) T 


Tf 1 (r^) = 


rl x V(r f ) T 

—i - (H rf H er V } 


(2.14a) 


(2.14b) 


The incident field is converted from spherical into rectangular 
components by using the following transformation equations: 


<H xf ’ H yf ’ H zf )T " D (H rf ’ H 0f’ V )T 


(2.15a) 


(E xf’ E yf ’ E zf )T = D (E rf* E 0f’ V )T 


(2.15b) 


Where D is defined by: 


D = 


sin cos <J>~ cos 0^ cos -sin 4y 

sin 0^ sin cos 0^ sin (jy cos <jy 
cos 0 r -sin 0 r 0 


(2.16) 


— O i 1 1 

Finally, the incident field ^ and in the feed coordinate system 
(Xf,yf,Zf) is transformed into the main coordinate system (x,y,z) by 
using the following equations: 


(2.17a) 


X 
LU 
1 

=T 

r E *fi 

E y 

= B 

V 

[ E z. 


L E zfJ 

PV 

=T 

r*fi 

H y 

= B 

V 

w 


[ H zfJ 


(2.17b) 
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Where B is given by 2.10. 

2.4 Aperture Integration Method 

The Aperture Integration Method is illustrated in Fig. 4. The 
surface is a Huygen surface enclosing the reflector and feed. 

This surface can be of arbitrary shape (typically spherical, cylindrical, 
etc). For the special case of a planar surface ^ ap (aperture plane) 
located in the radiation side of the reflector (the field over the 
shadow portion Ssh are assumed to be zero), the tangential fields 
on the aperture plane and the radiated far-fields are related by a 
Fourier transform (Ref. 6). The equivalence principle (Ref. 9) 
provides the basis for obtaining the radiated far-fields from the 
tangential 1 1 and it at ]£^ ap . Note that the aperture plane 
is taken perpendicular to the z-axi$. 

Let the tangential electric and magnetic field at X^ap be 
denoted by lf ap and lf a p respectively. Consider the following 
vector functions: 


'ap 


?<u,v> - t (x.y) e^ K<ux + v *>dx dy 


ap 


t(u.v) = // if (x,y) e^ K(ux + vy) dx dy 

v »P 


(2.18a) 


(2.18b) 


ap 


where u = sin 0 cos 4> 
v = sin 9 sin $ 


K = 2ir/X 




SCATTERED FIELDS 

e^d.Fu) 

OBSERVATION POINT 
(r,0,<J» 


4 Aperture and observation parameters for constructing the 
aperture Integration <AI) formulation 
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(r,0,4>) spherical coordinates of the far-field observation 
point. 

Since the aperture fields are tangent to the x-y plane, then 

^(u.v) = f x (u,v)>? + fy(u,v)y^ (2.19a) 

'j(u.v) = g x (u,v)lT + gy(u,v)”y (2.19b) 

The radiated far-fields from a current distribution can be obtained 
using potential theory (Ref. 10). In the aperture plane £»ap the 
induced currents are given by the following 

= n x H r (Electric current) (2.20a) 

x E r (Magnetic current) (2.20b) 

Where 1i r , ll r are reflected electric and magenetic fields (section 
2.6). The vector potential corresponding to the two induced currents 
(far-field approximation) are obtained from the following equations. 

t( r) = // (n x lt r )e - ^ K(r * r ’ } dx dy (2.21a) 

E ap 

?(r) = ^ f // -(n x lf r )e~^ K<r dx dy (2.21b) 

E ap 

where, and F* are the magnetic and electric vector potential 
respectively. 

r = x sin 9 cos <j> + y‘sin9sin4> + z N cos 0 (2.21c) 

’f , =x'T+yy' ( 2 . 2 1 d ) 

In the far-field region, 1: S and are related as in a spherical 
TEM wave. 


E s (r) 



( 2 . 22 ) 
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-v 

The far-fields E s (r) and H s (r) expressed in terms of vector potential 
functions /&r) and F*(r) are given by, 

^ S (r) = -jco &r) + juZ 'r xl^(r) 


^ S (r) = -jco F*(r) - ^-r x t( r) 


(2.23a) 
(2.23b) 

Equations (2.23a) and (b) imply that both electric (3*$) and 
magnetic currents (J* ms ) are included in the evaluation of far fields. If 
only the electric current or the magnetic current is present, the far-field 
electric field E* s (r) is determined by ^(r) or l*(r) alone. From 
the field equivalence principle (Love's equivalence principles) 
formulations, three different expressions arise for lf s (r), namely, 

( 1 ) Using both 


i s and J* ms ; 




E e = j K ! b7 - <f v cos 4) + f w sin 4) + Z cos 0 (g >( cos 4> - 9 V sin 4>)) 


4trr 


■iKe 


-jKr 


4ttc 


(2.24a) 


(cos 0 (f cos 4> - f sin 4>) - Z (g sin 4> + g cos 4>)) 
y x y x 


(2.24b) 


(2) Using 2J* S (Assume the aperture surface is a perfect 
electric conducting surface) 


r s _ iKZ e ^ Kr cos 0 
t 0 “ 


2-rrr 


,jKZ e~ ^ Kr cos 0 


2irr 


<9 y cos 4> - 9 X sin 4>) 


(g y sin 4> - g x cos 4>) 


(2.24c) 


(2.24d) 




(3) Using 2J ms (Assume the aperture surface is a perfect 
electric conducting surface) 


L 0 “ 2irr 


(f cos 4> + f sin 4>) 

x y 


(2.24e) 
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, il<p-j Kr cos 9 

E * " ^2^? (f y C0S ♦ - f x Sin ♦> (2 ' 24f) 

The results given by (2.24a) and (b) are the average of (2.24c) 
to (f) respectively. It requires the values of both tangential electric 
and magnetic fields at Sap- In practice it is more convenient to use 
either (2.24c) to (f) since they require the knowledge of only either 
the tangential magnetic or electric field over the aperture plane. 

This method is exact if the fields everywhere in Sq are known 
accurately. For electrically large reflectors the fields in the 
aperture plane will be usually small outside. an area defined by the 
projection of the reflector boundary on the aperture plane. 

2.5 Aperture Integration Method and the Fast Fourier Transform (FFT) 


To employ an FFT algorithm, the integrals in (2.18a) and (b) 
must be rearranged so that the form of integral matches the definition 
given by the FFT subroutines. In general, two dimensional FFT 
subroutines assumes the function, bei ng considered to be periodic in 
x and in y with period of 1 in both dimensions. Hence, a typical 
aperture plane grid as shown in Fig. 5, must be scaled accordingly. 

To illustrate the transformation of the given integrals into 
integrals with an FFT format consider the equation (2.18a) and extend 
the result to (2.18b). Each of these equations have two components, 

( f x > fy ) for Ed- (2.18a) and (g x ,gy) for (2.18b). First consider the 
component f x in Eq. (2.18a). 
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where E 3V is the x-component of E (x,y) in (2.18a). Using the 

d X dp 

following change in variable, 


x - xmin 
xmax - xmin 


(2.26a) 


y - ymin 
ymax - ymin 


(2.26b) 


and (xmax - xmin) = L x , (ymax - ymin) = Ly, Eq. (2.25) becomes. 


f „(u,v> 


where 


'1 PI 

, o Jc 


E (xmin + X (L ), ymin + Y (L )) 
3ix x y 


jK(u( L )X+v( L )Y) ) 
x e x y dX dY 


(2.27a) 


C 1 = (L x )(L y )e jK(u xmin+v ymin) (2.27b) 

The expression E fix (xmin + (l_ x )X, ymin + (L y )Y) within the 
integral may be interpreted as E (x,y) scaled in FE- -(a,b) within 

d X aX 

the limits 0 < a < 1 and 0 < b < 1 (Fig. 5b). FE ax (a,b) may be 
approximated by 


N2 M2 

FE ax (a,b) = ^ Y_ T mn e^ 2ir<mX+nY) (2.28) 

n=Nl m=Ml 


where T are the Fourier coefficients obtained by using an FFT 
mn 

subroutine. For example if one considers a 32 by 32 FFT grid points, 
then 


N1 = Ml = -15 and N2 = M2 = 16 
From (2.27a) and (b) and (2.28), 


f x (u,v) = C 1 


n 


0 J 


N2 

Z 

0 n=Nl 


M2 

Z T EXP (m.n.u.v) dX dY 
m=Ml mn 


(2.28a) 


where 
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jK((mX+u(L ) )X+(nX+v(L ) )Y) 

EXP (m.n.u.v) = e x y (2.28b) 

Interchanging the summation and integral signs and noting that 


,1 


IKcx 


dx = 


(s) 

P0 


(2.29) 


the expression for f is given by, 


f x (u,v> = C 1 


N2 M2 

2 c, 


. (m+n) 


n=Nl m=Ml 
where SF (m.n.u.v) is given by. 


mn 


SF(m,n,u,v) 


(2.30a) 


si 


SF (m.n.u.v) = 




(mX + ud^) 


si 


in (f~)(nX + v(L y ) )j 


(2.30b) 


(^(mX + u ( L x ) ) x(^-)(nX + v(L y )) 

In summary, to evaluate the integrals (2.18a) and (b) one 
needs to first calculate T mn of Eq. (2.28) by an FFT algorithm and 
then calculate f x (u.v) via (2.30a) and (b). The other components 
(fy,g x ,gy) Ec l s - (2.18a) and (b) can be easily transformed into an 
FFT format (Eq.( 2.30a) and (b)) following a similar approach. 

2.6 Computation of the Aperture Plane Tangential Fields 

The geometry for calculating the tangential fields on the 
aperture plane is depicted in Fig. 6. For a given feed point Pi and 
an observation point P 2 , a reflection point 0 r may exist on the 
reflector surface. This point is called a specular point (Ref. 34). 
This type of reflection satisfies Snell's law of reflection. The 
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6 Geometrical optics ray tracing for calculating the 
aperture fields 
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reflector surface is assumed to be a perfect conductor. Appendix B 
presents a method for obtaining the specular points on the reflector 
surface. The geometrical optics approximation of the reflected electric 
fields (Refs. 35 - 36) at an observation point P 2 is given by, 

t r <P 2 ) = Df e~^ d 2 (2(n • t 1 ) n -'th (2.31) 

where , 

d 2 distance between P 2 and 0 r 

incident electric field (see Section 2.3) 
n the reflector surface unit normal at 0 r 
DF divergence factor given by Eq. (2.32) 

DF = 1 ■ (2.32) 

x/l + (d 2 /Rj) \l 1 + (d 2 /R r 2 ) 

Rp R^ Principal radii of curvature of the reflected wavefront 
passing through 0 r . 

Appendix C contains a derivation of the principal radii of 
curvature of the reflected wavefront for a general surface of 
revolution. In summary for a given P] and P 2 , there may be more 
than one reflection point. Then the total reflected field at the 
aperture grid is the superposition of the contribution from each 
reflection point. If for any P 2 there is no specular reflection 
point on the the reflector surface, then the reflected field is zero. 

The approximation in Eq. (2.31) only takes into account the reflected 
rays. To improve the computation of the tangential fields on the 
aperture plane one may include edge diffracted rays (Refs. 37 to 43). 
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A similar expression to (2.31) can be found (Ref. 23) for this type of 
fields. In general edge diffracted fields and reflected field are 
superimposed at the aperture plane. For electrically large reflector 
antennas the contribution due to edge diffracted fields is very small 
and generally can be neglected. 

2.7 The Physical Optics Method 

The currents which excite the scattered field are induced on 
the reflector surface (assuming a perfect conductor) by an incident 
electromagnetic field (see section 2.3). Figure 7 illustrates 

the geometry for calculating the scattered fields using the physical 
optics approach. In the physical optics method the induced surface 
current distribution (Refs. 5 to 6) is approximated by 

? s = 2Cn‘ x ft’) On the illuminated side (2.33a) 

J* s = 0 On the shadowed side /L*sh- (2.33b) 

These approximations to the induced surface current are valid 
when the reflector size, the radii of curvature of the reflector 
surface, and the radius of curvature of the incident wavefront are all 
much larger than the wavelength. A fringing current component can be 
added to improve the current accuracy (Ref. 44). The fringing 
component is needed for applications that require the scattered fields 
to be known over wide observation angles away from the antenna 
boresight. 



REFLECTOR 



SCATTERED FIJELDS 

I s cn ,17(7) 

OBSERVATION POINT 
P: (r,e,<D) 


7 Geometry and parameters for constructing the physical 
optics (PO) formulation 
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The electromagnetic field solution to Maxwell's equations due 
to an electric current distribution in an unbounded region (Ref. 6) is 
given by Eqs. (2.34a) and (b). 

♦ K 2 f $ > | ^ dS (2.34a) 


E s (r) 


4ircoe 


*(r) = 


4-irojv. 


(jco^ x 1) | jKR dS 


(2.34b) 


^s 

Equations (2.34a) and (b) are valid at any arbitrary 
observation point P:(r,6,4>). In the integrands of Eqs. (2.34a) and 
(b) the operator sj operates on the coordinates of the source current 
alone. Expanding all the terms involved in the above integrands one 
obtains the following expressions. 



where 


(2.35a) 


(2.35b) 


K = — = u\f\Ic wavenumber 

R The distance from source element to the observation point P. 
A far-field approximation to Eqs. (2.34a) and (b) can be obtained by 
consideri ng 


R = r - "r • $ (for phase terms) (2.36a) 

R = r (for amplitude terms) (2.36b) 
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and by taking into account only 1/R variations in Eq. (2.35a) and 
(b). With these approximations, Eqs. (2.34a) and (b) are simplified 
into the following forms. 

d - (/ • r)r) e^ " r dS (2.37a) 
s s 

li s (r> . f| (^)' /2 <? s « r>e J ^* ~ r dS <2 ' 37b) 

where r is such that the following constraint is satisfied. 

r > — (2.38) 

D the diameter of reflector antenna 
A the operating wavelength 

The surface integration appearing in Eqs. (2.37a) and (b) is 
performed on the reflector curved surface. The differential area is 

(2.38) 

The square root factor above is commonly known as the Jacobian of the 
surface. It can be shown that the integral (2.37) and dS given in 
(2.38) may be represented in terms of a series of many two-dimensional 
Fourier transforms (Ref. 45). Typically, one needs the first few terms 
in this series to achieve a converging solution. The physical optics 
radiation integrals may be evaluated in many different ways (Refs. 12 
to 14) and its efficient evaluation has been a challenging problem. 


dS - *»\ j ' *( f ) 2 *(^) 2 


t 


<r> . 

4irr 
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2.8 Directivity Computation 

The far-zone electric field is usually divided into two 
orthogonal polarizations. Following Ludwig's definition 3 (Ref. 46), 
the following unit polarization vectors are introduced: 

R = 0 (a p e^P cos <j> + b p sin $) + $(-a p e^P sin $ + b p cos <j>) (2.39a) 
C = 0 (a p e~^P sin <J> - b p cos 4>) + $(-a p e"^P cos 4> + b p sin <j>) 

(2.39b) 

where (a p ,b p ,tp) are the far-field polarization parameters. These are 
related to the feed polarization parameters (section 2.3). For a 
single reflector these relationships are given by; 

a p — a, bp — b, ^ p = i|/ + tt (2.40) 

For example, a feed linearly polarized in the yf direction 
will produce a secondary radiation beam linearly polarized in the y' 
direction. A RHCP feed produces an LHCP secondary beam. For dual 
reflectors this relation may be more complex depending upon the 
coordinate systems and offset direction used. It is best to determine 
this relation in each individual problem.. 

If the secondary pattern can be expressed as (Eqs. (2.23a) and 
(b) and (2.37a) and (b) , 

* S< r) - ^ (e + * E +) (2 ' 4,) 

the reference polarization expression for ^ s (r) is given by (2.42a), 

ER(r) = • (R*") (2.42a) 

and the cross polarization expression is given by (2.42b) 

-4 

= 7 s • (C*) 


EC(r) 


(2.42b) 
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The directivity associated with the reference polarization is defined 
by 

DR (6, <j>) = ~~ (E s • R*) 2 r 2 /Z (2.43a) 

rad 

Similarly the directivity for the cross polarization is defined by 

DC (0, <|>) = p^- (^ S • C*) 2 r 2 /Z (2.43b) 

rad 

where P ra( j is the power radiated by the feed source (section 2.3). 

2 .-9 Numerical Results and Discussion 

This section presents some numerical results to establish 
numerical accuracy of the aperture Integration method (GO and GO + GTD) 
and the physical optics method. Far-field radiation patterns are 
presented for the following offset parabolic reflector configuration. 

Reflector geometry: D = 257.89 X diameter 

(refer to Fig. 8) f = 637.48 X focal length 

H = 135.51 X offset height 

Other parameters: Frequency = 30 GHz 

Aperture points = 200x200 
y'-polarized feed 
-18 dB edge taper 

A block diagram of a computer implementation (Refs. 22 to 23) 
for the aperture integration method (GO and GO V GTD) and the physical 
optics method are presented in Figs. 9(a) and (b) respectively. The 





32 



<a) Aperture integration algorithm. 

9 Block diagram of a computer implementation for calculating 
the secondary radiation pattern of a reflector antenna 





(b) Physical optics algorithm. 


Fig. 9. - Concl . 
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block diagrams describes all the computational routines that are 
required by both methods. 

Figures 10(a) and (b) show far-field radiation patterns 
(reference and cross polarizations) for the boresight case and 
Figs. 10(c) and (d) show a 6-beamwidth scanning case. As expected 
radiation patterns calculated by each method are in good agreement. 

The gain, half-power beamwidth (HPBW) and first sidelobe level for the 
boresight case are shown in Table 3. 


TABLE 3, - SUMMARY OF ANTENNA PERFORMANCE 


Method 

AI (GO) 

AI (GO + GTD) 

P0 

Gain (dB) 

56.88 

56.89 

56.85 

HPBW (deg) 

0.279 

0.281 

0.283 

SLL (dB) 

-33.89 

-35.2 

-33.3 


Note that the gain and HPBW obtained in all the methods are 
almost identical. The differences comes in the sidelobe level. The 
addition of the edge diffracted rays into the AI method improves the 
sidelobe level prediction, although for such electrically large 
reflector it does not seem to change the gain or the HPBW. 

Table 4 shows the computation times for the methods on an IBM 
370 computer. 


TABLE 4. - SUMMARY OF COMPUTER EXECUTION 

TIMES 


Method 

AI (GO) 

AI (GO + GTD) 

P0 

IBM 370. 

80 

150 

100 

C.P.U. 




(sec) 
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The AI (GO) method was the fastest as expected. The FFT 
algorithm makes this method numerically efficient. If the edge 
diffracted rays option is used then the AI become the slowest, but the 
most accurate. The PO is a trade-off between sidelobe accuracy and 
computational time. 
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CHAPTER 3 

REFLECTOR ANTENNA SURFACE ERROR REPRESENTATION 
3.1 Introduction 

Reflector antennas In a space environment are subject to 
continuous variation In temperature distribution, and are thus distorted 
from Its true geometrical shape (typically parabolic, hyperbolic, 
elliptical, etc.). The distorted reflector surface has in general a 
very complicated shape and hence can not be represented with an exact 
analytical expreslon. The analysis of a distorted reflector antenna 
defined by a set of discrete surface points require the use of numerical 
techniques. Many numerical techniques for analyzing the performance 
of reflector antennas defined by a set of points have been extensively 
reported In the open literature (Refs. 4, 32, 47 to 53). The most 
popular of these techniques represents the reflector surface either 
globally or locally by using polynomial splines (Refs. 54 and 55). The 
polynomial splines and other techniques require the surface points to 
be orderly labeled. This Is an undesired characteristic because it can 
result In a nonunique approximation for the desired reflector surface. 
The order of the spline polynomial necessary to best fit the reflector 
surface points Is in general unknown a priori, and therefore It 

Involves a trial and error procedure for checking the accurracy of the 

/ ' 

Interpolation. 
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The distorted reflector surface points can be best approximated 
by two analytical components, an undistorted surface component and a 
surface error component. The undistorted surface component is a best 
fit paraboloid polynomial for the given set of points and the surface 
error component Is the deviation of the actual surface points from the 
best fit paraboloid. This residual error component is then described 
with a sinusoidal Fourier series expansion. This approximation 
technique Is Insensitive to the labeling of the reflector surface 
points and can describe the surface errors very accurratly. Similar 
to a time signal the spatial spectrum of the surface error component Is 
unique to the reflector under consideration. Therefore spatial spectra 
can be utilized as a performance index for comparing distortion profile 
In reflector antennas. 

3.2 Description of the Problem 

The best analytical representation of the distorted reflector 
antenna surface points that uniquely Identifies surface errors can be 
obtained as follows: the reflector surface points are separated into 

two components, a best fit paraboloid component and a sinusoidal 
Fourier series expansion of the residual. Figures 11 illustrates a 
conceptual layout of the problem under consideration. 

In analyzing large reflector antenna performance it is 
necessary to accurately characterize the reflector surface points. 

Any derivation from Its Ideal geometry causes the antenna performance 
to degrade. The surface error component provides an independent 
performance index against which distorted reflector antennas can be 
compared. The surface error component by definition is a sinusoidal 
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11 Distorted reflector surface points separated into an undistorted 
surface component and a surface error component 
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Fourier series expansion of the difference between the actual reflector 

surface points and the best fit paraboloid geometry. Informations 

such as root mean square value (RMS), peak surface error and 

two-dimensional distortion profile can be calculated from the surface 

error component. The Fourier coefficients In the series expansion 

represents the spatial spectrum that uniquely Identifies the distorted 

reflector under study. The best fit parabolic surface represents the 

reflector antenna surface In an average sense. When the surface error 

Is zero the best fit paraboloid surface reduces to the Ideal or the 

design surface geometry. 

3.3 Undistorted Surface Component 

A generalized two dimensional second order polynomial 

(paraboloid of revolution) Is given by Eq. (3.1). 

2 2 2 

S (x ,y,z) = a 1 x + a 2 y + a 3 z + a^ xy + a g yz + a g zx + a 7 x 

+ a g y + z + a ]0 (3.1) 

This polynomial form can be transformed Into a different coordinate 

system In which Eq. (3.1) will have the following form, 

S' (x ' ,y ' , z 1 ) - x' 2 + y’ 2 - FL z' (3.2) 

Where FL Is a constant. Figure 12 shows a generalized 

parabola of revolution (Eq. (3.1)), the (x,y,z> coordinate system and 

the (x'.y'.z 1 ) coordinate system. These two coordinate system are 

related by set a of Eulerian rotations (Ref. 56) and a translation. 

Let the distorted reflector surface be represented by 

OS (x^.y^.z^) = 0, 1=1 , N (3.3) 

Using the method of least squares (Ref. 57) the distorted 

reflector surface points from an approximating parabolic polynomial 

given by Eq. (3.1) can be determined in the manner described below. 
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The minimization index is defined as 


I = £ (DS(x. ,y. ,z. ) - S (x.,y.,z.)) 2 
i=l ill 

where 

DS (x.,y.,z.) nth order surface polynomial representing 
distorted reflector surface points. 

S (x^.y^.Zj) the approximating parabol ic polynomial in 
Eq. (3.1) 

N total number surface points 

Substituting Eq. (3.1) and (3.3) into (3.4) yields the following 
result. 


r.t 


1 = 1 


(a, x* * a 2 yf ♦ a 3 z? * a„ x,y, * a 5 y,z, ♦ a 6 z,x, 


+ a 7 x i + a g + z i + a io )< 
For minimum I it is required that, 


9 1 3 1 3 1 


3 I 


3 a 1 ' 3 a 2 3 a 3 


3 a 


= 0 


10 


Equation (3.6) results in the following set of equation. 


3 I 
3 a. 


N 


Z 2 <S(x. ,y. ,Zj )>x^ ■ 0 


3 I 
3 a. 


= £ 2 (S(x.,y.,z.))y 2 = 0 


i = l 


3 I 
3 a. 


£ 2 (S(x.,y z,))z 2 = 0 
i=l liii 


3 I 
3 a. 


£ 2 (S(x. ,y. ,z.))x i y. = 0 


(3.4) 

the 


(3.5) 

(3.6) 

(3.7a) 

(3.7b) 

(3.7c) 


(3 . 7d) 
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8 I 
8 a. 


- £ 2 (S(x 1 ,y ) ,z,»y f z ) - 0 


U1 


<3.7e) 


8 I 


3 a 


= 21 2 (S(x . ,y. ,z.))x.z. = 0 


6 i = l 


(3. 7f ) 


H;- £ 2 <S(x | .y ) .z,))x, -0 


(3 . 7g) 


3 I 


8 a 


= £ 2 (S(x. ,y. ,z. ))y. = 0 


8 I 


8 i = 


N 


A 2 < S(X 1’ y i’ Z i> )1 =o 


(3. 7h) 


(3.71 ) 


l 10 i =1 

Equations (3.7a) to (i) can be arranged in a matrix form given by, 


R (50 = E 


(3.8) 


where R is a 9x9 real symmetric matrix with the following 


coeff i 

cients 

X . 

y.x. 

X 

N 

Vi 

y . z . x . 

z . x . 

x . 
l 

y.x. 

x , 
1 


1 


1 1 

iii 

i i 

l l 


x.y. 

y. 

2-y: 

z.y. 

X . Z V 

x . y . z . 

x.y. 

y. 

y. 


i i 

1 

1 1 

1 1 

1 1 1 

i i l 

i i 

1 

i 


X . 2 . 

y 

Z . 

x . y . z . 

y.z. 

x . z . 

x . z . 

y.z. 

z . 


1 1 

i i 

1 

i i i 

1 1 

l i 

1 1 

i i 

1 

N 

x.y. 

y.z. 

z . x . y . 

x.y. 

x . y . z . 

x . y . z . 

x.y. 

x.y. 

x.y. 

R =Z 

i i 

i i 

i i i 

1 1 

l i i 

ii i 

i i 

i i 

i i 

i = l 

x . y . z 

y.z. 

y . Z ; 

x . y . z . 

y.z. 

x . y . z . 

x . y . z . 

y.z. 

y.z. 


v i i 

i i 

i i 

i i i 

i i 

ii i 

l i l 

i i 

i i 


X . z . 

x . y . z . 

x . z . 

x . y . z . 

x . y . z . 

x . z . 

x . z . 

x. y. z. 

X . z . 


1 1 

i i l 

i i 

ii i 

l l l 

i i 

i i 

ii i 

1 1 


X . 

x.y. 

x . z . 

x.y. 

x.y . z . 

x . z . 

x . 

y.z. 

X . 


1 

1 1 

i i 

i i 

l i i 

i i 

i 

i i 

1 


x.y. 

y. 

y i z ; 

x.y 

y.z. 

x . y . z . 

x • y • 

y. 

y. 


1 1 

1 

i i 

i i 

i i 

i i i 

i i 

i 

i 


X . 

y. 

z . 

x.y. 

y • z • 

x . z . 

x . 

y. 

i 


1 

i 

i 

i i 

l l 

i i 

l 

i 



(3.9) 
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E is a column vector with the following elements. 


E = 


i = l 


Vi 


Vi 


2 i 


Vi z i 


V 


x i z 


x i z 


V 


z i 


(3.10) 


and X column vector with the following elements, 

'*1 


X . 


L a ioj 


(3.11) 


The solution to Eq. (3.8) provides the optimum coefficients for 
a best' fit paraboloid defined by Eq. (3.1) in the least square sense. 

Differential geometry (Refs. 58 and 59) forms the basis for 
obtaining the translational and rotational relationship between the 
(x,y,z) coordinate system and the (x',y',z') coordinate system. In 
order to find the translation between both systems it is necessary to 
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R - a] x 2 + a 2 y 2 + 


find the vertex of the paraboloid defined by Eq. (3.1). Let's define 
the following variables, 

P = a 3 (3.12a) 

0 = as y + ag x + 1 (3.12b) 

a 4 xy + 37 x + ag y + a^g (3.12c) 
and rewrite Eq. (3.1) as a function of the new variables P,Q and R. 
This yield the following, 

S ( x , y , 2 ) = Pz 2 + Qz + R (3.13) 

the roots of the Eq. (3.13) describes a generalized paraboloid in the 
form z = f(x ,y) , 


dSL 


i 




4PR) 


2P 


(3.14) 


the positive root is discarted because yields a solution in the 
negative z-direction. Differential geometry formulae can be used on 
Eq. (3.14) for obtaining the vertex of the paraboloid. The vertex is 
located at a point (xO.yO.zO) at which the two principal curvatures 
have equal magnitudes. The principal curvatures of a surface of 
revolution are defined as follows: 

K i = K m + v /(K m 2 - K g ) (3.15a) 

K 2 = K m -\/<K m 2 - K g ) (3.15b) 


where 


K = Eg- 2f F + eG ( 


mean curvature) 


2( EG - F c ) 


K_ = 


r ^ 

eg - f 


(Gaussian curvature) 


(3.15c) 


(3.1 5d) 


(EF - F ) 
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In Eq. (3.15c) to (d) the variables (e,g,f) and (E,G,F) are 
the first and the second fundamental coefficients associated with a 
surface of revolution. These coefficients are defined as follows. 


E = 1 


(3.16a) 

F 

' f ‘ f y 

(3.16b) 

G = 1 


(3.16c) 

e 

■ 4f xx 

(3.1 6d) 

f 

■ Af xy 

(3. 16e) 

g 

' 4f yy 

/ j \ 

(3.1 6f ) 

A = - 


(3. 1 6g) 

f 

_ df (x,y) 

(3.1 6h) 

x 

d x 


f xx 

d 2 f(x,y) 

(3.161) 

d 


f y 

_ d f(x,y) 
d y 

(3. 16j) 

f 

d 2 f(x,y) 

(3.16k) 

yy 

dy 2 


f = 

d 2 f(x,y) 

(3.161) 

xy 

d x d y 



with these definitions, the required derivatives are obtained from 
Eq. (3.14) with following results, 

f _ _ I_ U2 _ ]_ d V(Q 2 - 4PR) 
x ~ 2P d x 2P dx 


(3.17a) 
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f 1 d 2 0 1 d l J (0 2 - 4PR) 

dx d x 

(3.17b) 


f 1 d 0 1 d\/(Q 2 - 4PR) 

r y " " 2P d y " 2P dy 

(3.17c) 


f 1 d 2 0 Id 1 J (0 2 - 4PR) 

f yy = - 2P d y 2 - 2P d y 2 

(3. 17d) 


f 1 d 2 0 ' 1 d 2 y (Q 2 - 4PR) 

r xy 2P d x dy ~ 2P d x d y 

(3. 17e) 

where 



<Q 2 - 

4PR) = x 2 (a 6 2 - 4 a 3 a|) + y 2 (a 5 2 - 4 a 3 a $ ) 



+ xy (2 a^ a^ - 4 a^ a^) + x (2 ag a g - 4 a^ 

a 7 > 


* y <2 a 5 a, - ( a, a g ) , U g 2 - 4 a 3 a, 0 > 

(3.18) 

Using 

the following substitutions 



m, - <a 6 2 - 4 a 3 a,) 

(3.19a) 


”2 - <a 5 2 - 4 a 3 a 2> 

(3.19b) 


m 3 - (2 a 5 a 6 - 4 a 3 a„> 

(3.19c) 


= (2 ag a g - 4 a 3 a^) 

(3. 19d) 


m 5 • <2 a 5 a, - 4 a 3 a 8 > 

(3. 19e) 


m 6 ‘ <a 9 2 - 4 a 3 a 10 > 

(3. 19f ) 


the required derivatives in Eq. (3.17a) to (e) are obtained as follows, 

T x = d<9 Vx 4 PR) = n» 1 x h- y m 3 ♦ w» 4 (3.20a) 

T d 2 (Q 2 - 4 PR) 


<3. 20b) 
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r _ d <Q - 4 

y ~ d y 


PR) 


= y + x m^ + m^ 


t d 2 (0 2 - 4 PR) 

1 , = = = m„ 


yy 


d y 


T d 2 (0 2 - 4PR) _ m 

T xy d x d y m 3 


(3.20c) 


(3.20d) 


(3.20e) 


By letting V = n /cq 2 - 4PR) , one obtains the following result. 


O.0.5V-" 2 T x 


(3.21a) 


^2 = 0.5 V 1/2 T + -0.25 T 2 V" 3/2 
d x x 

0 = o.5V 1/2 T 
d y y 

^2 = 0.5 V~ 1/2 T + -0.25 T 2 V' 3/2 

d y^ yy y 

rrdV °- 5 ' r ' /2 T xy + - 0 ' 25 T X T y v " 3/2 

H - a6 


(3.21b) 
(3.21c) 
( 3 . 2 1 d ) 

(3. 21 e) 
(3.21f) 


0,-0 

d x^ 

d_fl _ a5 

d y - a5 


(3. 21 g) 


(3 . 2 1 h) 


= 

d / 


= 0 


(3.211) 


d -Q _ 0 

d x d y 


(3.21 j) 


Equations (3.21a) to (j) evaluates the set of derivatives required in 
Eq. (3.14) for calculating the vertex. The location of the vertex is 
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found numerically. Reflector surface points are searched for a unique 
point at which the two principal curvartures are equal (K] = < 2 >. 

The rotational characteristics between the (x,y,z) coordinate 
system and (x'.y'.z 1 ) coordinate system can be obtained by first 
finding the inward unit normal at the vertex point. This will be the 
z' axis direction. The other two axes directions relative to the 
(x,y,z) coordinate system are found from the normal direction. The 
unit normal at the vertex is given by. 


f ,» * V + V 

<f 2 . f 2 , f 2 > 
x y z 


= z , 


(3.22) 


where f x and f are defined in (3.16h) and 3 . 1 6 j ) respectively, and 
ri f ( x v x ) 

f = — • The projection of this unit normal n on the x-y 

plane is given by, 


xy 


V * V 




(3.23) 


The Eulerian rotations can be calculated by using the following 
expressions. 


cos (Gl) --n • <-y> (3.24a) 

cos (G2) = n • z (3.24b) 

G3 = 0 

The surface points defined in the (x,y,z) coordinate system are 
transformed into the (x'.y'.z 1 ) coordinate system using the following 
matrix equation: 
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e 11 e 1 2 e l 3 X 1 ” *0 
e 2 1 e 22 e 23 y i " y 0 
e 31 e 32 e 33. Z i " Z 0 

where 


ei i = cos 9 

(3.26a) 

e ^2 = sin <j> 

(3.26b) 

ei 3 = 0 

(3.26c) 

e 2 1 = - cos 9 sin <j> 

(3.26d) 

e 22 = cos 6 cos <t> 

(3.26e) 

e 23 = sin 6 

(3.26f) 

e 3 i = sin <J> sin 0 

■(3.2 6 g ) 

632 = - sin 6 cos 4 > 

(3.2 6 h ) 

633 = cos 6 

(3.26i) 


With ail reflector surface points transformed to the 
(x'.y'.z 1 ) coordinate system the best fit paraboloid of revolution 
have the desired form given by Eq. (3.2). The only parameter unkown 
in Eq. (3.2) is the constant FL. 

Let's define the following least square index 

N P 
I'= £. (DS’(x' y' xl) - S'(x' y' z\)r (3.27a) 

i=l iii iii 

where 

DS' (x‘ ,y' ,z' ) nth order surface polynomial representing the 
distorted reflector surface points. 

S' (x'-,y' ,z' ) paraboloid of revolution defined by Eq. (3.2) 

substituting Eq. (3.2) into (3.27a) we obtain. 



( i=l ,N) (3.25) 
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.•.I 

, i = l v ■ 


- z i. 


Differentiating with respect to 1 / FL one obtains, 

dl' 


N 


= 1 


0 = 2 
i = l 


d 1 /FL 


(xl) 2 + (y ! ) 2 


FL 


= 0 


- Z i 


((xj) 2 , <yj> 2 ) 


Solving for 1 / FL yields 


1 /FL = 


iii z i( (x ; >2 - ( yj >2 ) 

N 2 2 

ZL ( (x ! ; + . (y! ) > 

i = 1 1 1 


Thus the equation for the paraboloid of revolution in the 
(x'.y'.z 1 ) coordinate system is given by 

x 1 2 + y' 2 = ( 1 / F L ) z ' 

3.4 Surface Error Component 

The surface deviations of the distorted reflector surface 
points with respect to the best fit paraboloid is defined as the 
surface error component. The surface error component can be 
expressed as, 

zz: = Dz! - z: , i = 1 ,N 

where 

( x l ,y ! ,Dzi ) reflector surface points 

2 2 

x! + y! best fit paraboloid 

7 • _ — ! — — 

z i 4F 


1.27b) 


(3.28) 


(3.29) 


(3.30) 


(3.31) 


(3.32) 
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A two-dimensional continuous function like the one above and 

j 

defined in the domain, -K < x < K and -L < y < L can be approximated 
by a two-dimensional Fourier series (Ref. 60). 



b cos(-) si n(^y^) 
mn K L 


„ . , mirx . „ / 1 

sin(-^-) cos(- 


sinC—) 


sin(^) 


(3.33) 

where the Fourier coefficients a mn , b mn , c mn and d mn are given by, 

'L PK 

a mn = f(x ’ y) cos cos (^) dx dy (3.34a) 

-L -K 

'L fK 

b mn = f<x,y) cos(^) cos(~ y ) dx dy (3.34b) 

-L -K 

'L fK 

c = f(x,y) cos( n] 7 ^) cosC^) dx dy (3.34c) 

mn J K L J 

-L -K 
'L PK 

d mn = f(x,y) cos(^) cos(~ y ) dx dy (3.34d) 

ll 

If the function z = f(x,y) have the following symmetry 
conditions, 

f(-x,y) = -f(x,y) and f(x,-y) = -f(x,y) (3.35) 

Then Eq. (3.33) can be simplified to, 

00 CO 

f(x,y) = 2. £ d mn sin(^) s i n( n f^) (3.36) 

m=l n=l m ^ L 


a two-dimensional sinusoidal Fourier series. 
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In general the surface error component zz’ is defined over a 

i 

region xmin < x'< xmax and ymin < y'< ymax (Fig. 13). Equation 

(3.36) requires zz' to satisfy the symmetry conditions given by Eq. 

i 

(3.35). A coordinate system (x",y",z") is used for rearranging the 

reflector surface points in a format such that the symmetry conditions 

are satisfied. First the surface error points in Eq. (3.2) are 

translated to the first quadrant of the (x",y",z") coordinate system. 

The other three quadrants are fil led with data points correspondi ng -to 

the symmetry conditions in Eq. (3.35). This process is illustrated in 

Fig. 14. The Fourier coefficients of the zz" surface error 

i 

component can be obtained as follows. 


mn 


4 

K L 


A A 

x y 


£ 


1=1 


M 

£ 

j=l 


zz 


,mir 


ntr 


( x ,yV) s\nC~r‘) ■sin(~-y'j) (3.37) 


where 



M x - 7- ♦ ' 

(3.38a) 


X 




(3.38b) 

A x’ A y 

sampling spacings in x and 
y dimensions respectively 

(3.38c) 

xf " x i 

- xmin translation in x dimension 

(3. 38d) 

yv = y; 

- ymin translation in y dimension 

(3.38e) 


Let us define the Fourier matrix D, 



(XMIN,yMAX) 


(XMAX,yMAX) 
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13 Reflector surface points defined over a rectangular 


grid 




14 Illustration of reflector data points augmentation and 
. translation 
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\ 


11 

d !2 

d l 3 d ln 

21 

d 22 

••• d 2n 


(3.39) 



Where d mn are the Fourier coefficients defined by Eq. (3.37). 
The elements of the matrix D defines the spatial spectra of the 

surface error for a given reflector antenna surface. The order n of 

\ 

this matrix is estimated by including only the coeffients with 
amplitude greater than X/200 or it can be arbitrarily chosen. 

Another parameter of interest is the root mean square (RMS) value of 
the surface error component, which is defined as, 

RMS = (d 2 + d 2 + ... d, 2 + dJ + d 2 + ... + d 2 
II I i. In 21 22 2n 

22 2 1/2 
+ ... + d n f + d n , + ... + d ^ Y ,C (3.40) 
nl nZ nn 

3.5 Numerical Results and Discussion 

Figure 15 represents a block diagram of a computer simulation 
of the above stated problem. The reflector surface points are usually 
obtained from a holographic, photogrametri c or any other surface 
detection technique. The computer algorithm was tested with known 
distortion profiles and TRASYS-SYNDA-NASTRAN (Ref. 4) simulated 
thermal distortions superimposed into a reflector antenna geometry. 
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The distortion profiles considered for the simulation are 
described in Table 5 and the reflector geometry illustrated in Fig. 16. 
Case A in Table 5 describes a small (X/20 or less) distortion profile 
and case B illustrates a large (of the order of several wavelengths) 
distortion profile. The frequency used for the simulation was 10 
GHz. A set of equally spaced data points (100x100) were needed to 
generate the distorted reflector surface. 


TA8LC 5. - OISTORHON PROFILE PARAMETERS 



Fourier Coefficients 
(In wavelength) 

<>ii 

d, 2 

d 13 

"21 

*22 

<*23 

d 31 

<*32 

d 33 

Case a 

♦0.0500 

♦0.0700 

-0.0300 

♦0.1500 

-0.0040 

♦0.0900 

-0.0033 

-0.0083 

-0.0310 

Case B 

♦0.0300 

.♦0.3000 

-2.000 

-3.000 

♦0.4000 

♦3.000 

-1 .000 

♦5.000 

-0.3000 


The estimated surface errors and best fit parabolic surface geometry are 
described in Table 6. 


TABLE 6. - ESTIMATED DISTORTION PROFILE PARAMETERS 



Fourier Coefficients 
(In wavelength) 


d 12 

<*I3 

d 21 

<*22 

d 23 . 

d 31 

<*32 

<*33 

Case A 

♦0.0497 

♦0.0710 

-0.0305 

♦0.1501 

-0.0039 

♦0.0913 

-0.0032 

-0.0082 

-0.0312 

Case 8 

♦0.0318 

♦0.2999 

-1.998 

-2.999 

♦0.4108 

♦2.988 

-0.9999 

♦5.100 

-0.2991 


Estimated focal length: 
Case A 3.290109 ft 
Case B 3.313091 ft 


Figures 17(a) and (b) show the estimated surface profile 
corresponding to case A and case B and their respective far field 
radiation patterns. The results presented in Table 6 are in good 
agreements with the input distortion by Fourier coefficients on 
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0.000250 
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<a) Case A (small distortion profile). 
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E-plane radiation pattern and their estimated distortion 
profile 
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Table 5. It was found that the Fourier matrix had higher order 
coefficients with non zero values but they were one order of magnitude 
less than the lowest amplitude of the Fourier coefficient presented on 
Table 6. 

The technique was compared with results obtained by using a 
spline polynomial technique for approximating the distorted reflector 
surface points (Ref. 4). In this case the distorted reflector surface 
points were obtained by simulating thermal deformation with a 
TRASYS-SINDA-NASTRAN computer programs. In brief, TRASYS and SINDA 
are used to characterized the on-orbit thermal envlroment, NASTRAN 
calculates the thermally Induced mechanical distortions. Figure 18 
shows the temperature distribution on the reflector antenna surface 
for the case under consideration. The reflector geometry Input to the 
thermal programs Is presented In Fig. 19. The frequency considered 
was 28.75 GHz and a 100x100 surface points were used for the analysis. 
The far field radiation pattern corresponding to the thermal 
slmulatatlon case Is presented In Fig. 20. The continuous line pattern 
corresponds to the polynomial spline algorithm (Ref. 4) and the dotted 
line pattern the corresponds to the best fit paraboloid and Fourier 
series expansion approximation. There Is no major differences between 
beam direction or sldelobe levels. Indicating a good agreement between 
the two techniques. The best fit paraboloid and Fourier series 
algorithm was very slow, about 3 hr of c.p.u. time In an IBM 370 
computer. The spline polynomial algorithm takes about 1/2 hr on CRAY 
xmp computer. The long computation can be justified as a trade off to 
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obtain vital Information about the distorted surface characteristics. 
These are the amplitude spectra of the surface error, the distorted 
surface profile, the RMS value and largest deviation on the reflector 
antenna surface. These are not directly available by any of the other 
existing techniques. 

3.6 Concluding Remarks 

One advantage of the developed technique is that it can be 
easily Implemented In any existing reflector antenna secondary pattern 
computational methods. It can easily be extended to nonparabolic 
reflectors surfaces (spherical, planar, hyperbolic, elliptical, etc) by 
modifying the least square polynomial approximation. In applications 
Involving the fabrication and design of precision reflector antennas 
the technique can be used as a computer aided tool. Information such 
as the average focal length, a root mean square of the surface error 
(RMS), surface error profile and the amplitude spectra for the reflector 
antenna under consideration are easily accessible from the algorithm. 

One draw back, of the algorithm Is that It Is slow. This can be 
Improved by using faster computer such as a CRAY and optimizing the 
computer codes for faster computing time. 
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CHAPTER 4 

COMPENSATION METHOD I: TRANSMITTING APPROACH 

4.1 Introduction 

The performance characteristics (gain, sidelobe level, etc.) 
of space communication reflector antenna systems degrades as the 
reflector surface distorts due mainly to thermal (Ref. 4) effects 
caused by solar radiation. The compensating technique described here 
will maintain the desired boreslght directivity (antenna performance) 
Independently of thermal effects on the reflector surface. A feed 
array can be designed to compensate for thermally induced surface 
distortions. 

For a better understanding of reflector surface error 
compensation, it is convenient to think of the antenna in the receiving 
mode. Figure 21(a) shows a uniform plane wave Incident on a perfect 
parabolic reflector. In the geometrical optics sense, all rays will 
converge to a single point at the focus of the parabola, where they can 
be collected by a single feed element. Figure 21(b) shows a uniform 
plane wave incident on a distorted parabolic reflector. In this case, 
the rays will not converge to a single point, but will spread over the 
focal plane. Reflector surface error compensation will be achieved, 
if a feed array can capture all the energy spread over the focal plane 
and sum it constructively in the feed network. Improvement in the 
antenna performance due to compensation in the receive mode will 




21 Plane wave incident on a reflector antenna 
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(b) Rays spread over the focal plane of a distorted reflector surface. 
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produce Improvement In the antenna performance In the transmitting 
mode by reciprocity. 

The concept of conjugate field match (CFM) is utilized to 
determine the complex excitation coefficients of the feed array. 
Specifically the compensating feed array excitation coefficients are 
determined from the complex conjugate of the focal plane electric 
field (Ref. 33) distribution due to a uniform plane wave incident on 
the distorted reflector surface from the desired beam direction. In 
the transmitting approach (Ref. 61) the focal plane electric field 
distribution is determined by computing the far zone electric field in 
the desired beam direction due to a dipole source at a specified 
location in the focal plane. This computation procedure is performed 
for all specified feed array element locations in the focal plane and 
their conjugate values are used as the compensating feed array 
excitations . 

4.2 Determination Of The Compensating Feed Array Geometry 

It is assumed that the distorted reflector surface (usually 
described by a discrete set of points) is known a priori. The 
distorted reflector surface points can be approximated with two 
analytical surface components (see chapter 3): an undistorted surface 

component and a surface error component. The undistorted surface 
component is mathematically expressed with a best fit paraboloid 
polynomial and the surface error component is described as a 
sinusoidal Fourier series expansion of the residual. The surface error 
component with the largest spatial frequency and amplitude product 
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determines the feed array area required (Ref. 62) to capture ail the 
energy spread in the focal plane (Fig. 21(b)). 

The geometry of the best fit paraboloid (undistorted surface 
component) is illustrated in Fig. 22. This geometry is completely 
specified by the following angles: 

©off : off-set angle of the reflector 

e max : half-angular extent of the reflector surface 

The minimum area consisting of D by D (Fig. 22) in the 
focal plane needed for capturing about 80 percent of the energy spread 
(Ref. 62) is approximately given by 

D * 1T5-F- <’ * cos W <4 ' ,) 

ma x 

where e = 2.71828. . . , 

B: amplitude of the largest surface error component 

T: number of periods of the largest surface error 

component over the extent of the reflector surface. 

The required number of elements that can fit into the 
compensating area D (Eq. (4.1)) varies according to the desired array 
configuration and on the element aperture area. Typically horns, 
microstrips and dipole antennas are considered as radiating elements 
in either a hexagonal or rectangular array geometries. The gain and 
sidelobe level requirements of the subject reflector antenna will 
usually determine the number of radiating elements on the feed array 
(Refs. 63 and 64). A rule of thumb criteria for the element spacing 
for a given feed array configuration is to have minimum spacing in 
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22 Best fit parabolic geometry 
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order to reduce the spillover when illuminating the reflector. This 
is a good criteria for cases In which the feed array is located in the 
near field zone of the distorted reflector. For cases in which the 
feed array is located in the far field zone, the spacing should be 
determined such that there are no grating lobes generated within the 
reflector surface. 

4.3 Computation of the Compensating Feed Array Excitation Coefficients 

The focal plane electric field distribution contains the 
necessary information for finding the required number of elements, 
their best location and their proper complex excitations. Conceptually 
this field distribution results from an incident uniform plane wave 
impinging (Fig. 23) on the distorted reflector surface from the desired 
beam direction. The complex element excitation are obtained by taking 
the complex conjugate of the focal plane electric field distribution. 
This concept is known as conjugate field matching principle (CFM). 

Let's assume that the feed array contains N elements located 
in the focal plane. The complex excitation coefficients, represented 
by a column vector I, 

I = til, I 2 , I 3 W 1 • <4.2) 

are to be determined. 

Let us first calculate the far field radiation pattern of the 
distorted reflector antenna due to the mth element in the feed array 
(Fig. 24), namely due to the following array excitations, 
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23 Uniform plane wave illuminating the distorted reflector 
antenna 
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element located in the focal plane 
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I m = 1 , and all other 1} = 0 (4.3) 

The corresponding far zone electric field can be written as, 

— J K r ^ ^ 

tjr) = - (RE m (e,<|)) + CF m ( 6 , 4 >)) (4.4) 

m r m T m 

where (r - ®, 9 = 4> = ^BEAM^ s P herica1 coordinate of the 

a -a 

desired observation point in the far field zone, (R,C) are unitary 
complex vectors describing the reference and the cross polarization of 
the reflector antenna (Ref. 46) respectively. 

Let the column vector, 

^^BEAM’^BEAM^ = ^1 (4.5) 

represent the reference polarization electric field in the observation 

direction OgEAM’^BEAM^ clue each the individual elements in the 
feed array. The far fields in Eq. (4.5) can be easily obtained by 
using either a physical optics or an aperture Integration technique 
(see chapter 2). Typically the physical optics method is chosen due 
to its computational ease. This method requires that the induced 
current distribution on the distorted reflector surface be integrated 
in order to obtained the radiated electric field at the desired 
observation direction ■ ^BEAM ^ ' 

Appendix D shows that by applying Lorentz reciprocity 
principle (Ref. 9) the focal plane electric field distribution of a 
distorted reflector antenna can be indirectly obtained from the 
transmitting electric fields due to dipole sources located in the focal 
plane (Refs. 61 and 62). In the receiving solution the focal plane 
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electric field distribution Is directly obtained by integrating the 
induced current distribution due to an Incident uniform plane wave on 
the distorted reflector. 

In the transmitting approach, the compensating feed array 
excitation coefficients are determined by using the reciprocity 
principle. The resulting complex excitation coefficients are given by 
Eq. (4.6). 

1 = E * (e BEAM ,< ^BEAM > (4 ‘ 6 

where ^gEAM’^BEAM^ the des,red beam direction. This approach for 

determining the compensating array excitation coefficients is called 

indirect conjugate field matching (ICFM). This method has been widely 

used for more than 10 years (Refs. 11, 33 and 66 to 77) for 

applications requiring compensation of degraded antenna performance 

due to off-focus feeds. 

4.4 Numerical Results and Discussion 

A block diagram of a numerical implementation for the 
transmitting approach for calculating the compensating feed array 
excitations is presented In Fig. 25. Notice that for each of the feed 
element location (x m ,y m ,z m ) in the focal plane a far zone electric 
field is computed (Eq. (4.4)) in the desired beam direction 
^BEAM’^BEAM^ using a physical optics algorithm. 
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25 A computer implementation of the ICFM 
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As an example of a surface error compensation using a feed 
array,, consider the undistorted reflector and a given surface error 
profile presented in Figs. 26(a) and (b> respectively. The operating 
frequency is 30 GHz and the feed element type are considered to be 
square horns with a (1.6X x 1.6X) aperture dimension. The undistorted 
E and H-plane far field reflector antenna patterns with the single 
feed horn located at the focal point are shown in Fig. 27. The 
undistorted antenna boresight directivity is 58 dB. With the given 
surface error profile (Fig. 26(b)) superimposed on the undistorted 
reflector geometry (Fig. 26(a)) the resulting E and H-plane far 
field reflector antenna patterns with a single feed horn located on the 
focal point are shown in Fig. 28. The antenna boresight directivity 
is seen to be reduced. to 38 dB (20 dB loss). It is required to 
compensate for the loss in directivity due to the surface distortions 
with a feed array. The first step is to design a compensating array 
configuration that will capture most of the energy spread in the focal 
plane. By applying Eq. (4.1) to the reflector geometry and considering 
the distortion profile shown in Figs. 26(a) and (b), we obtain the 
following compensating area, 


D = 


(e)(ir)(0.7697)(0. 5)X 
sin (20°) 


(1 + cos (43°)) = 17. 3X 


The radiating elements are square horns with a (1.6X x 1.6X) 
aperture dimension. The number of elements necessary for a rectangular 
array geometry are estimated to be 121 elements (a 11 x 11 array) at 
about 1.6X spacing. The corresponding compensated E and H-plane 
reflector antenna pattern are shown in Fig. 29. The compensated 





T = 0.7697 
B = O.S\ 


(b) Surface error profile 


26 Distorted offset reflector antenna geometry 
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boresight directivity is within 0.969 dB of the undistorted boresight 
directivity; clearly showing a good design for the feed array size and 
element excitation coefficients. For comparison purposes let us consider 
another feed array geometry shown in Fig. 30. Note that the majority 
of the array elements are located along the vertical dimension of the 
array and very few elements are allocated horizontally. This feed 
array design has been made by taking advantage of the a priori 
knowledge of the distortion profile. Since there is no distortion 
variation in the horizontal dimension of the reflector it is expected 
that elements located horizontally in the feed array will contribute 
very little to performance compensation. The corresponding compensated 
E and H-plane antenna patterns are shown in Fig. 31. The compensated 
directivity for this case is within 1.7 dB of the undistorted boresight 
directivity. No major differences in the compensated E or H-plane 
pattern shape or boresight directivity is observed when compared to 
those obtained by using a 11 x 11 rectangular array configuration. The 
reason for this results is due mainly to the one-dimensional variation 
of the surface error profile (Fig. 26(b)). Therefore the excess of 
radiating elements in the 11 x 11 array distributed horizontally in 
the feed array did not add any extra compensation to the antenna 
pattern shape and directivity. 

In general the surface error profiles due to thermal effects 
are estimated to have two dimensional variations (Ref. 4). This 
implies a large number of radiating elements to be required for a 
rectangular feed array configuration. Consequently the distortion 
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error profile Is a factor that needs to be taken into consideration 
when designing the compensating array. This technique clearly shows 
an improvement to the degraded reflector antenna directivity. The 
only draw back is lack of sidelobe control in the compensated antenna 
performance. For applications involving space reflector antenna where 
the main concern lies in maintaining the boresight directivity 
independently of thermal effect the transmitting approach certainly 
lend itself to an adaptive implementation. 
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CHAPTER 5 

COMPENSATING METHOD II: RECEIVING APPROACH 

5. 1 Introduction 

Array-fed reflector configurations are very desirable for 
space communication multiple and contour beam antenna systems. The 
radiation performance of these reflector antennas is degraded because 
of surface distortions induced by thermal effects from solar radiation. 
Distortions caused by thermal gradients are inherently a large scale 
phenomenon. By using computer programs such as TRASYS-SYNDA and 
NASTRAN (Ref. 4) to simulate thermal distortions, antenna performance 
degradations can be calculated. For example, boresight directivity 
loss of 3 to 8 dB, sidelobe level increase of 10 dB, and pointing loss 
of 0.1° to 2°, have been predicted by TRASYS-SINDA and NASTRAN 
simulation. The chapter describes a technique to compensate such a 
degradation, so that the antenna performance is virtually independent 
of thermal effects on the reflector surface. 

The compensating technique is briefly .outl ined below. Normally 
in a multiple beam antenna system, each beam is excited by a single 
element in the feed array (one-to-one excitation) such an excitation 
is not good enough if the reflector surface distortion is present. In 
the latter case, we will use a cluster of N elements (7 to 35 or even 
more feeds) to excite a single beam (N-to-one excitation). By adjusting 
the complex excitation coefficients of the cluster feed, the antenna 
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degradation due to reflector surface distortion can be partially 
restored. The degree of restoration depends on the number of feeds In 
the cluster, and the precision of realizing the desired excitation. 

Obviously, the N-to-one excitation scheme Is more complex than 
the conventional one-to-one scheme. The former requires both the 
excitation amplitude and phase control of all elements in the feed 
array. With the advent of monolithic microwave integrated circuits 
(MMIC) , such a control becomes realizable without excessive effort 
(Ref. 78). It should be mentioned that in addition to electromagnetic 
compensation techniques studied here, other efforts have been carried 
out to reduce the impact of thermal distortion. Most of them are of 
mechanical nature (Refs. 79 and 80). 

5.2 Problem Description 

Given the desired boresight directivity, far field sidelobe 
envelope and distorted reflector surface points (x,y,z), the 
compensation problem is to determine the minimum number of feed array 
elements, their best location and their excitations that will give rise 
to the desired antenna performance despite of the surface distortion. 

The approach here is based on the receiving antenna mode, as 
outlined below. Let the distorted antenna be illuminated by an incident 
"plane wave" from a prescribed direction and with a prescribed 
polarization. That direction and polarization are those of the main 
beam when the antenna is in a transmitting mode. The "plane wave" has a 
planar phase front as in the ordinary plane wave case, but with a 
tapered amplitude, in contrast to the constant amplitude in the ordinary 
plane wave case (Fig. 32). 




32 Distorted reflector antenna being illuminated by a tapered 
plane wave 
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Using this concept, the resulting fields In the focal plane of the 
reflector antenna can be calculated by Integrating the Induced current 
distribution on the distorted reflector surface. The focal plane 
electric field distribution contains the necessary Information for 
determining the required number of feed array elements, their best 
location and complex excitations necessary to compensate for the 
surface errors. The set of excitation coefficients of the 
compensating feed array Is obtained as a complex conjugate of the 
received focal plane electric field distribution. With this set of 
excitation coefficients the feed array Illuminates the distorted 
reflector surface with a wavefront designed to compensate for 
reflector surface distortions. This method of obtaining the 
compensating array excitations is called direct conjugate field 
matching (DCFM). Next, the difference of the present DCFM from the 
other compensating technique will be explained. 

In the past, conjugate field matching principle was used for 
obtaining the compensating element excitation coefficients in 
applications where degraded reflector antenna performance arises from 
off-focus feeds (Refs. 11,76,77, and 82). Recently the same technique 
was applied to the compensation of large scale surface distortions 
(Refs. 61,65, to 68) with a *feed array. It is named indirect 
conjugate field matching (ICFM) because the transmitting mode and 
reciprocity theorem were used In determining the feed array 
excitations. The difference between DCFM and ICFM lies in the 
fol lowl ng facts : 
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(i) DCFM is based on the receiving antenna mode. Hence the 
tapered plane wave can be conveniently added in the 
formulation. Since the aperture field and the far field 
of an antenna are related by a Fourier transform, there 
is an explicit relation between the plane wave taper and 
the antenna sidelobe level. 

(ii) The compensating excitation coefficeints in ICFM are 
invariably determined under the condition that the 
reflector antenna is in a transmitting mode. The 
aperture field taper is not an explicit design 
parameter, and therefore cannot be conveniently 
controlled. 

As a consequence, we note that ICFM often leads to antenna 
patterns with high sidelobe envelopes (see Chapter 4). Such a problem 
can be alleviated by the present DCFM. 

5.3 Computation of the Focal Plane Fields 

Consider a linearly polarized (y> tapered plane wave incident 
on the distorted reflector surface from an (Fig. 33) observation 
direction (0,4)). 

t 1 ( x . ,y 1 , z 1 ) = E(p) e + ^ Kr * * u (5.1) 

Where , 

E(p) = c + (1 - c> C 1 - (p/a) 2 ] (5.2a) 

p - \/[(x s } 2 ♦ (y s - y c ) 2 ] <5. 2b) 
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c = 10 ET/2 °; ET: is the edge taper of 

the incident plane wave 

y c ,a: reference lengths (Fig. 34) 

( V y s ,z s > : distorted reflector point 

(5.2c) 


^ >*v 

r = x $ x + y s y + z s z 

(5. 2d) 


K = 2 ir/X: wavenumber 

(5.2e) 


X: operating wavelength 


II 

sin © cos 4> x + sin 0 sin <J> y + cos © z 

(5.2f ) 

<0 , 4>) : spherical angles describing the 

observation direction. 


/■s 

U = 

A A 

u X + u V + u 7 z: polarization vector 

x y z 

(5.2g) 


2 

u = ( - sin © cos 4> sin <t>) 
d 

(5.2h) 


2 2 2 
u = (cos © + sin © cos <i>) 

^ d 

( 5 . 2 i ) 


u_ = (-cos 0 sin © sin <j>) 
z ■ d 

(5.2j) 

2 

{(sin © cos 4> 

sin 4>) 2 + (cos 2 © + sin 2 © cos 2 4>) 2 



2 1 /?. 

+ (cos 0 sin 0 sin 4> )} 

(5.2k) 


One needs to calculate the electric field at a point 
(xi,yi,Zj) in the focal plane of the distorted reflector antenna. 
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The incident electromagnetic field E\ induces an 
electric surface current distribution on the distorted reflector 
surface, which is assumed to be a perfect conductor. The induced 
electric current distribution (Ref. 5) is given by 

3 - 2 n x ft 1 (5.3) 

where 



surface normal 


(5.4) 


z = f ( x , v )• distorted reflector surface 
s s’-'s 

By assuming also that the incident field is a transverse 
electromagnetic wave (TEM) one can make use of the following 
auxi 1 iary relationship, 


where 

Z = 120 ir: free space wave impedance 

s | = - sin 9 cos 4> x — sin 0 sin 4> y - cos 9 z (5.6) 

Substituting (5.6) into (5.1), yields 

3 S = 2 {n x [s. x t 1 ]} (5.7) 

The scattered E and ft fields from the distorted reflector surface 
in terms of the electric current surface distribution J* s is given 
(Ref. 6) by 
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t . j-=i- 

4iro> e 


{K 2 3 s «|» + (3 $ • $)$t> dS 


(5.8a) 


r r 



$t> dS 


<5. 8b) 


Where 


Co = Z-af: and f is : the operating frequency, in Hz 

-12 

e = 8.854x10 :j permitivity of free space (farads/meter) 



(5.9a) 


r . \/{(x i - x $ ) 2 + (y. - y $ ) 2 + (z. - z,.) 2 } (5.9b) 

The time dependent factor has been taken as e3“t and is 
omitted through out the analysis. 

In the integrands of Eqs. (5.8a) and (b), the operator 1 
acts locally on the the coordinates of the source element. By 
expanding Integrand of Eqs. (5.8a) and (b) we obtained the following 
expressions. 


% («|r> 


• (* • 


(5.10a) 


(3 • $) $ (»| ») 


t3 s • V'l + r * ?) A * 


-r (*♦?). 




(5.10b) 


where r^ = 


(x 1 - x $ ) x + (y t - y $ ) y + (z. - z $ ) z 


(5.10c) 
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Substituting (5.10a) and (b) into Eqs. (5.8a) and (b) gives the 
corresponding focal plane fields E^P and H^P 


f P 


<x i '^i * 2 i ) " W 


K 2 3 - K 2 
s 


( 3 s * ~ r \) 


r l * 


J. 


? ♦ *) ( j s • MJ 4- 

^‘VXi-V-fefJ [(V'i) (jk + 1) ] ^ dS 


dS 


(5.11a) 

(5.11b) 


5 . 4 Determination of the Compensating Feed Array Geometry and 
Element Complex Excitation Coefficients 

This section will address the determination of the complex 
excitation coefficients of the compensating feed array . But before 
the analysis is presented a broad guideline will be established for 
obtaining the best compensating array configuration. Several 
procedures for designing an array-fed reflector antennas have been 
reported in the literature (Refs. 23,33,61,63,64,68). The most 
popular approach utilizes the focal plane field distribution as a 
guideline for obtaining the array boundaries. 

To provide reflector distortion compensation, a feed array must 
be designed such that it will capture the energy spread in the focal 
plane (receiving mode) and sum it constructively in the feed network. 
Improvements in antenna performance due to compensation in the receive 
mode will produce improvements in antenna performance in the transmit 
mode by reciprocity. To determine the necessary compensating feed 
array area for capturing the energy spread in the focal plane fields 
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it Is necessary to truncate the focal plane fields to within a 
predetermined relative value. This is necessary because in general 
this fields cover a large focal plane area. Another aspect that should 
be considered in estimating the compensating feed array area lies in 
the electronics of the devices of the feed network. For example the 
limiting factors in the amplifiers of the feed network are mostly due 
to finite dynamic range, noise level and finite number bit quantization 
levels. These limiting factors determine the minimum amplitude level 
that can be resolved by the feed network. The minimum amplitude level 
is used to locate the compensating feed array boundaries by filtering 
out focal plane fields lower than the minimum amplitude. This process 
geometrically maps the boundaries of the compensating feed array. 

With boundaries of the compensating array estimated by truncating the 
focal plane fields, the next step is to choose a feed array element. 

The element type is usually determined by the desired reflector gain 
and sidelobe level (Ref. 63). The number of elements in the feed 
array can be estimated from the truncated focal plane area divided by 
the aperture area of the individual element. The element spacing 
should be keep to a minimum to minimize the spill-over losses. Once 
the compensating feed array configuration is determined, the relative 
amplitude and phase of each array element must be determined. 

Let us consider a distorted parabolic reflector, antenna with a 
feed array of N elements. The complex excitation coefficients, . 
represented by a column vector I, 
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I = CIi . 12. 1 3 , - - - 1 N ] T <5.12) 

are to be determined. Let us first calculate the focal plane electric 
field distribution (5.11a) at each (xj.yj.Zj) feed array element 
location. The focal plane electric field distribution is denoted by, 

£ fp <x ? ,y i , z 1 ) = Ex[ P x + EY^ P y + EX[ P z 1=1, N (5.13) 

Let's define the unit vector p^ to be the focal plane 
reference polarization. The distorted reflector geometry and the 
vector p f are presented in Fig. 35. The unit vector is given 
by Eq. (5.14). This vector is directed from the focal point to the 
geometrical center of the distorted reflector. 

V x + y r y + <z r ■ F) 2 

S = — (5.14) 

vK 2 ♦ > c 2 * <z c - f>2) 

where 

(x ,y z )• coordinate of the geometrical center of the 
distorted reflector surface 

F: focal point 

By considering only y polarized feed elements, the corresponding 
reference polarization unit vector is given by, 

P f = P x y + P z z (5.15) 

where 

-x 

p = — c — (5.16a) 

V “c 2 * ♦ <z c - F)2> 

(ZC - F) 

x ‘ > - nh 


(5.16b) 
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35 Geometry for calculating the focal plane reference polarization 
unit vector 
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Considering only the focal plane electric field in the direction of the 
focal plane reference polarization vector p^., 

E. = ? fp (x.,y.z.) • ; f (5.17) 

Let the column vector 

E = [E lt E 2 ,E 3 E n ] T (5.18) 

represent the focal plane reference polarization of a distorted 
reflector at each element location in the feed array. 

In the transmitting mode the distorted reflector antenna is 
being illuminated by the compensating feed array with excitation 
coefficients that are obtained as a complex conjugate of the reference 
polarization electric field (Eq. (5.18)) in the focal plane. The 
resulting compensating complex excitations in the feed array are given 
by 


■r 


L U 

>2 


E 2 

‘3 


E 3 



, e n. 


Equation (5.19) describes the direct conjugate field match 
(DCFM) approach or the receiving approach for compensating the large 
scale surface distortions. It is called direct because it computes 
the excitation coefficients directly by integrating the induced 
currents on the distorted reflector surface. 
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5.5 Numerical Implementation 

A block diagram of a numerical implementation of the receiving 
approach for compensating large scale surface distortions is presented 
in Fig. 36. It consists of three numerical algorithms; a surface 
analyzer, a compensation analyzer and a radiation pattern analyzer. 
These algorithms form a complete set of techniques for compensating 
and evaluating the large scale surface errors (thermal distortions 
type) in a distorted reflector antenna systems. 

The first algorithm, reflector surface analyzer is Illustrated 
in a block diagram in Fig. 37. It takes the distorted reflector 
surface points as the input and finds the best fit paraboloid 
polynomial (see chapter 3). The difference between the best fit 
paraboloid polynomial and the input points (residual) is expanded into 
a sinusoidal Fourier series. Essentially this algorithm converts the 
input discrete points into two analytical surface components. 

The second algorithm, the compensation analyzer is illustrated 
in a block diagram in Fig. 38. It determines the focal plane electric 
field distribution due to a tapered plane wave Incident on the 
distorted reflector from an observation direction. The algorithm uses 
physical optics techniques for obtaining the focal region fields. The 
compensating feed array are obtained by assigning the complex conjugate 
values of the focal region field at the corresponding feed element 
locations. In summary this algorithm will compute the relative complex 
excitation coefficients of the feed array. 
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36 A computer implementation of surface error analysis and 
compensation algorithm (receiving approach) 
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« 


37 A block diagram of the reflector surface analyzer algorithm 











38 A block diagram of the reflector surface compensation 
algorithm 
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The third algorithm, distorted reflector radiation analysis is 
illustrated in a block diagram in Fig. 39. It determines the 
compensated boresight directivity and far field radiation pattern. 

The array-fed distorted reflector geometry and compensating array 
excitation coefficients are the input to this algorithm. The algorithm 
utilizes a physical optics approach (chapter 2) for obtaining the far 
field radiation pattern and boresight directivity. 

The three algorithm are combined to produce the generalized 
compensation algorithm for large scale surface distortions. Basically 
two inputs are to be defined, the feed array geometry and the distorted 
reflector surface points. The feed array geometry usually contains the 
feed element type (typically microstrip, open waveguide, horn, etc.), 
array shape (typically rectangular, hexagonal, etc.), number of 
elements and relative location of each elements. The distorted 
reflector surface points are usually obtained from optical, microwave 
holography, photogrametry, or any other metrology technique. The 
distorted reflector surface points can be described in either an 
equally spaced grid or in a nonequally spaced grid. 

5.6 Numerical Results and Discussion 

To evaluate the usefulness of the receiving approach for 
compensating large scale surface distortions with a feed array, the 
distorted reflector surface will be simulated by superimposing a given 
distortion profile (Fig. 40) into an offset parabolic reflector 
(Fig. 41). First the effects of truncating the focal plane fields to 
a finite amplitude level on the compensated boresight directivity and 
sidelobe envelope will be investigated. Next the effects of tapering 
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39 A block diagram of. the distorted reflector secondary 
pattern algorithm 









110 

the Incident plane wave amplitude on the compensated antenna 
performance will be studied. For comparison purposes the result of 
the conjugate field matching technique applied to a simulated distorted 
reflector will also be presented In this section. 

The parameters describing the simulated surface error profile 
(Fig. 40) are outlined In Table 7. The surface error profile was 
chosen to resemble a typical large scale distortion profile caused by 
thermal effects. The undistorted and distorted E-plane radiation 
pattern obtained by using a single feed element are shown In Figs. 42 
and 43 respectively. Table 8 presents a description of the feed 
parameters. The calculated undistorted and distorted boreslght 
directivities are 39.1 dB and 30.30 dB respectively. The undistorted 
and distorted antenna performance provides the upper and lower limits 
In which the compensation technique can be applied. 

The effects of truncating the focal plane fields at a finite 
amplitude level are evaluated by first computing the focal plane 
electric field distribution due to a tapered plane wave (-17 dB) 
incident on the distorted reflector surface from the boreslght direction 
The -17 dB taper In the Incident plane wave is the required amount of 
taper for achieving a first sidelobe level of -33 dB in the compensated 
pattern. This value Is calculated from the aperture illumination edge 
taper consideration that will give rise to a far field sidelobe level 
of -33 dB. Figure 44 shows a three-dimensional plot of the focal 
plane plane electric field distribution resulting from Illuminating 
the distorted reflector with taper plane wave (-17 dB) from the 
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CONTOUR CONTOUR 
LABEL VALUE 

A 0.008000 
B 0. 004000 
C 0.002000 
D 0.001000 
E 0.000500 
F 0.000250 
G 0.000125 
H 0.000000 
I -0.000125 
J -0.000250 
K -0.000500 
L -0.001000 
M -0.002000 
N -0.004000 
0 -0.008000 


X, METERS 


40 Reflector surface error profile 





41 Undistorted offset parabolic reflector geometry 





TABLE 7. - SURFACE ERROR PROFILE PARAMETERS 


TWO-DIMENSIONAL FOURIER COEFFICIENTS IN WAVELENGTHS 

d n 

d 12 

d 13 

d 21 

d 22 

d 23 

d 31 

d 32 

d 33 

0.001500 

0.002100 

-0.000900 

O.OORSOO 

-0.00012R 

0.002578 

-0.000099 

-0.000250 

-0.000928 

TYPICAL ELEMENT IN THE SERIES: d mn SIN [(y - 0. 1)(m)(2rt] SIN [(x + 0.5)(n)(27T) 
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TABLE 8. - SUMMARY OF FEED SOURCE PARAMETERS 


OPERATING FREQUENCY, GHz 10 

FEED POLARIZATION y LINEAR 

E-PLANE FEED PATTERN C0S qE (0) . . . qE = 17 
H-PLANE FEED PATTERN COS qH (0) . . . qH = 17 
FEED LOCATION (FOCUS), METERS . . (0,0, 1.0) 
EDGE TAPER, dB -17 
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boresight direction. Three different focal fields cutoff (truncating) 
amplitude levels, -10, -20, and -30 dB will be evaluated. Figure 45 
shows a center cut plot of the focal region fields and also shows the 
corresponding array areas (D1 , D2, D3) for the given focal plane field 
cutoff levels. Figures 46(a) to (c) presents rectangular feed array 
geometries corresponding to the -10, -20 and -30 dB focal fields cutoff 
levels respectively. The element spacing and element pattern are the 
same for all cases. The parameters describing the array elements are 
presented in Table 9. The resulting compensated E-plane radiation 
pattern are shown in Figs. 47(a) to (c) respectively. The calculated 
compensated boresight directivity were 36.51, 34.7, and 32.3 dB 
respectively. Note the sidelobe envelope was best compensated for the 
-20, and -30 dB cases. The best compensated boresight directivity was 
2.59 dB from the undistorted reflector with a single feed case. 

In evaluating the effects of tapering on the incident plane 
wave the compensating rectangular array configuration that corresponds 
to the -30 dB focal fields cutoff level (Fig. (46(c)) will be selected. 
The taper levels that will be examined are 0, -9, -17, and -26 dB. For 
each taper amplitude level a new set of feed array complex excitation 
coefficients are computed. Figures 48(a) to (d) presents the 
compensated E-plane radiation patterns for the 0, -9, -17, and -20 dB 
taper level case respectively. The calculated boresight compensated 
directivi ty were 38.87, 37.56, 36.51, and 32.12 dB respectively. 

Notice that the best compensated boresight directivity occurs for the 
0 dB (Fig. 48(a)) taper, but this case also produced the worst 
compensated sidelobe level envelopes. The sidelobe envelope is best 
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44 Three-Dimensional plot of the focal plane electric field 
distribution 


RELATIVE AMPLITUDE, dB 
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45 Center cut plot of the focal plane electric field 
distribution 
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Fig. 46. - Concl . 


TABLE 9. - SUMMARY OF FEED ARRAY ELEMENT 
PARAMETERS 


OPERATING FREQUENCY, GHz 10 

FEED POLARIZATION y LINEAR 


E-PLANE FEED PATTERN COSTCO) . . . qE = 1 
H-PLANE FEED PATTERN COS qH (0) . . . qH = 1 
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47 Compensated E-plane radiation pattern 


RELATIVE AMPLITUDE, dB 
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Fig. 47. - Cont. 
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(c) -30 dB focal field cutoff level case. 


Fig. 47. - Concl . 
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compensated for an amplitude taper of -17 dB (Fig. 48(c)). For the 
case of -9 dB taper a high sidelobe envelope was obtained and for the 
-26 dB taper good sidelobe level was achieved but with a lower 
compensated boresight directivity. This result indicates that 
tapering the incident plane wave provides a direct control on the 
compensated sidelobe structure and also achieves a good boresight 
compensated directivity (wi thing 1 to 3 dB of the undistorted boresight 
directivity). 

In the previous analysis the receiving approach for 
compensating large scale surface errors has been demonstrated. The 
next step is to compare the receiving compensating technique with 
results obtained by the other conjugate field matching approach (see 
chapter 4). Let us consider the undistorted offset parabolic geometry 
and a surface error profile shown in Figs. 49 and 50 respectively. 

The parameters describing the simulated surface error profile are 
presented in Table 10. Figures 51 and 52 show the undistorted and 
distorted E-plane radiation patterns obtained with a feed element at 
the focal point. Table 11 presents the parameters describing the feed 
element. The antenna boresight directivity has been reduced to 38 dB 
with the distortion. Let us consider the compensating array geometry 
in Fig. 53. The compensated E-plane radiation pattern using the 
receiving approach is shown in Fig. 54. The taper level of the 
incident plane wave is -15 dB. The compensated E-plane radiation 
pattern using the transmitting approach (see chapter 4) is shown in 
Fig. 55. The transmitting approach assumes a uniform plane wave 



127 



4UI — — J- 1 1 1 i I i i I i 

-10 -8 -6 -4 -2 0 2 4 6 8 10 

6. DEG 


<a) 0 dB amplitude taper case. 


48 Compensated E-plane radiation pattern 
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e, DEG 

<c) -17 dB amplitude taper case. 


Fig. 48. - Cont. 
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Fig. 48. - Concl . 
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49 Offset parabolic reflector geometry 



AZQ) 
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T = 0.7697 
B = 0.5A 


50 Surface error profile 
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incident on the distorted reflector. The compensated directivity for 
both approaches were within 2 dB of the undistorted value. The major 
differences occurred in the compensated sidelobe level envelopes. Note 
the sidelobe level for the receiving approach were 3 to 5 dB lower than 
those produced by the transmitting approach. This again indicates that 
amplitude the tapering of the incident plane wave have a direct control 
on the compensated sidelobe envelopes. 

Considering the computational speed, the receiving approach was 
fast and numerically efficient in computing the compensating feed array 
coefficients. The receiving approach requires N steps for computing 
the distorted reflector currents induced by the tapered plane wave. The 
transmi tting approach wi 1 1 require M*N steps, to calculate di storted 
reflector currents Induced by the M feed elements in the focal plane. 

This makes the transmitting approach slower and numerically inefficient. 

For the case that has just been illustrated the receiving algorithm was 
very fast requiring about 20 min of C.P.U. time about 1 hr of C.P.U. 
time in an IBM 370. 

The receiving approach for compensating large scale reflector 
surface distortions was demonstrated and compared with other exciting 
techniques. The effects of truncating the focal plane fields to a finite 
amplitude level and the effects of the tapering the incident plane wave 
in the receiving approach were investigated. The results shows that the 
receiving approach was able to compensate sidelobe levels and also 
provided a reasonably good compensated directivity (withing 1-3 dB) . The 
method also indicates computational ease and suggests its application in 
an adaptive type of implementation involving large reflector antenna system. 
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52 Distorted E-plane radiation pattern 
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53 Compensating array geometry <37 elements) 
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TABLE U. - SUMMARY OF FEED ARRAY PARAMETERS 

OPERATING FREQUENCY, GHz 30 

FEED POLARIZATION .... LINEAR 

E-PLANE PATTERN COS qE (0). qE = 20 

H-PLANE PATTERN C0S qH (0). ...... qH = 20 

FEED LOCATION (FOCUS), meters . . . (0,0,3.35) 
EDGE TAPER, dB . . . . -11 
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54 Compensated E-plane radiation pattern using the 
receiving approach (DCFM) 
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CHAPTER 6 


SUMMARY 

6.1 Concluding Remarks 

The reflector antenna surface distortion arised from thermal 
effects can be adequately compensated by using an active phased array 
feed. There are two approaches to design the array excitation 
coefficients, the indirect conjugate field matching (ICFM) discussed 
in chapter 4 and the direct conjugate field matching (DCFM) in chapter 5. 
The DCFM method produces the lowest compensated sidelobe and the ICFM 
method produces the highest directivity. Due to its computational 
ease (computational speed) DCFM lends itself to a simple realization 
in hardware (if the distorted reflector surface is known at all time), 
and may be very useful in an adaptive type implementation. This is 
particularly true considering the recent advances in MMIC feed array 
technology, where each radiating element in the feed array can be 
controlled independently to provide the compensating excitations 
(ampl 1 tude and phase) . 

6.2 Suggestion for Future Work 

The compensation technique should be experimentally 
demonstrated by using a distorted reflector with an MMIC feed array. 

A detection scheme is necesary to obtain the distorted discrete surface 
points as input to the compensation algorithm for any real time 
application. The two algorithms can be used in an adaptive 
implementation for compensating thermal distortion in a space 
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environment. The method presented here can be extended for 
nonparabolic surfaces, such as spherical reflectors and should find 
applications in antennas used In satellite communications, space 
radiotelescopes, microwave power transmission and radiometers. 
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APPENDIX A 

POWER RADIATED BY AN ARRAY FEED 

The generalized feed array geometry Is Illustrated In Fig. A-l. 
The feed position corresponding to the mth element Is at (x m ,y m ,z m ) and 
Its complex current excitation Is I m . The radiated electric field 
due to this element in the array can be expressed by, 

e -JKr . 

E „< r > - — V9-*> <*•'> 

where the vector function f^(0,<{>) defines the mth element polarization 
and far-fleld radiation pattern. This vector function may be 
approximated by, 

frn(e,<t>) = 'eUE m (9)(aeJ'l' cos <J> + b sin <}>) + $UH m (0)(bcos <|> - ae^ sin 4>) 

(A. 2) 

where 

UE m (0) E-plane pattern of the mth element 
UH m (6) H-plane pattern of the mth element. 

Typically, the far-fleld element patterns may be approximated by a 
(cos 0) ^ function, i .e. , 

UE m (0) = (cos 0>Q E (A. 3a) 

UH m (0) = (cos 0>9 H (A. 3b) 

The radiated far-fleld of the feed array Is the product of the element 
pattern given by A.l and the array factor, i.e.. 




A_1 Geometry of the 


generalized feed 


array 
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hAi) - l v* 1K " * *' 

m=l 


(A. 4) 


where E m <r) is given by A.l and 

l? = sin 0 cos 4> + sin 0 sin 4> y + cos 0 “z (A. 5a) 

*p|n = x m x + y m y + z m z (A. 5b) 

The total time-averaged radiated power of the array is given by, 

>2* 


"rad Z Q 


it/2 

0 


r 2 stn9d9d* <A.6> 


where is the electric field radiated by the feed array given 

by A. 4 and Z = \fyTTc = 120 ir. 


Substituting (A. 4) into (A. 6) we obtain 


rad 


M M 

- L I 

m=l n=l 


'II* — 

\V» I n 

V/2 

1 0 J 

0 


INTEGRAND (m,n,0,(j>)) d0 d<fr' 


where 

jKu • (f -i ) 2 

INTEGRAND (m,n,0,<j>) = E m (r) • E*(r) e m r c sin 0 

Defining the power radiated as follows, 


(A. 7a) 


(A. 7b) 


rad 


M M 

2 VS V, 

m=l n=l 


(A. 8) 


where A is given by the expression in the brackets of Eq. (A. 7a). 
mn 

It can be shown that, 


?{( 


E m (r) • E* (r) = { (a t UE m (0)UE n (0) + b t UH m (0)UH*(0)j cos c <j> 


m 


+ (b 2 UE (0)UE*(0) + a 2 UH(0)UH 
m n m 


(■ 


(" 


; <e) ) 


m 


sin 2 4> 


+ ab cos <|) <UE m (0)UE*(0) - UH (0) - UH m (0)UH 
1 m n n m 


*(0)^ sin 2$J 


sin 2$ ) (A. 9) 
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substituting (A. 9) into (A. 7b), then the expression for k is 

mn 


given by 
1 


mn 2 


0 J 


fir/ 2 
0 


f2ir 

0 


a 2 UE (0)UE*(0) + b 2 UH m (e)UH*(0)}cos t * + 
m n m n 1 T 








+ lb 2 UE(0)UH*<0) + a 2 UH (0)llH*(0)] sin 2 * + ab cos u; (UE (0)UE*<0> 
1 m n m n ' T 1 m 


Y (l 


- UH m (0)UH*(0)) sin 2 4 .) x e 


m n 

jKr mn <s1n 6 - s1n © cos (<J> - <j, mn ) 


mm 


) 

+ cos 0 m „ COS 0) 
mn 


sin 0 d 0 d<j> 


where 


r = / ( x - x ) 2 + (y_ - y n ) 2 + (z - z m ) 2 

mn ^ m n n m m 


r m J n 


"m m 


8 mn = cos 
mn 


-1 


- z n) 


mn 

♦ = tan"' ^ 

¥ mn x m - x„ 

m n 

With the use of the following integration formulas. 


f2ir j 5 cos < 4> - <t>) 

e c d 4 > = 2 ir J Q (^> 


r mn 


(A. 10a) 


(A. 10b) 


(A. 10c) 


(A.IOd) 


(A. 11a) 


„2i r 


J O 


2.1 j E cos < 4 > - 4> r ) 



d<j) = ir 0 Q (£) ? ir cos (2 

(A. 11b) 


The result after <$> integration is reduced to 
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A - -2L 
mn ' Z 


f ir/2 


UE (0)UE*(0) + UH <0)UH*(0> J n <Kr sin 0 mn sin 0) 
m n m n / 0 mn mn 


jKr cos 0 m „ cos 0 

x sin 6 e 00 + ( 1 — Cq) J 


ir/2 


UE (0)UE*(0) 
m n 


- UH <©)UH*<0) J 0 (Kr sin 0 m „ sin 0) sin 0 e 
m n / 2 mn mn 


jKr cos 0__ cos 0 i* 
J mn mn d0 


where 


< 2— 72a) 


C n = 1 + 2ab cos u/ sin 2 <p „ + (a 2 - b 2 > cos 2 4> 
0 T mn 


r mn 


(2— 12b) 
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APPENDIX B 
SPECULAR REFLECTION 

For a given feed at point Pi (xi.yj.zi) and an observation 
point P 2 <x 2 .y 2 .z 2 ) <F1g. B.l), a reflection point 0 r may exit on the 
reflector surface <x,y,f(x,y)). The vectors 

d i = < x - X]) x + (y - y"i> y + (z - Z]) z (B.la) 

d 2 * <x - X 2 > x + <y - y 2 > y + (x - X 2 ) z (B.lb) 

are the connecting vectors between PI and 0 r , and from 0 r to P2, 
respectively. A necessary condition for which the reflection points 
must satisfy is that the total distance (dl + d2) must be stationary, 
i .e. , 

^ <d ] + d 2 ) = 0 and ^ (d ] + d 2 > = ° <B.2) 

after differentiation we obtain. 


1 

d l 

<x- 

' X 1 

) + (f(x,y) 

- 

V 

df 

dx 

1 

(x - 

X£> + (f(x,y) 

- z >£ 
z 2 ; dx 

= 0 












(B.3a) 

1 

d l 

<y- 

- *1 

) + (f(x,y) 

- 

V 

df 

dy 

1 

<y - 

y 2 > + (f(x,y) 

- z 
z 2 ay 

= 0 












(B.3b) 



A 

root (x.y.z 

= 

f ( X , 

.y)> 

of the 

two 

nonl inear equations i n 


(B 

,3a) 

and 

(b) gives 

location of 

a reflect! 

ion point. For 

a given 

point 

Pi 

and 

P2. 

there may 

be 

none, one, or 

more 

than one reflection point. 

It 

may 

be 

shown that 

(B. 

,3a) 

and 

(b) is 

equi 

ivalent to sati 

sfying Snell 's 


law of reflection (Ref. 6). 
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B.l Geometry for the specular reflection calculation 


r 
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The system of Eqs. (B.3a), and (b) can also be satisfied if 
Pi, 0 r , and are col linear Suich an unwanted root may be 
eliminated by imposing an additional condition. 



(B. 4) 

where e is a small positive number. One can preset e * 0.00001. A 
root of (B.l) may or may not fall inside the boundary of the reflector. 
Each root (x,y,f(x,y)) must be checked if is inside the desired 
boundary. 
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APPENDIX C 

PRINCIPAL RADII OF CURVATURES OF THE REFLECTED WAVEFRONT 


Differential Geometry (Ref. 83) forms the basis of the 
computation of principal curvatures of the reflected wavefront 
(RpR^) . The three orthonormal base vector of the incident ray bundle 
are chosen to be (Fig. C-l) 


~1 

1 


y x x. 


i 


y x x. 


x(z - z 1 ) - z(x - x 1 ) 

\/((z - z } ) 2 + (x - x^ 2 ) 

~i ~1 
x 3 * X 1 


(C. la) 


(C.lb) 


x(x - x.) + y(y - y.) + z(z - z,) 
x 3 - ■■ " ■■■ ' ■ ■ — (C.lc) 

\/(x - x^ 2 + (y - y^) 2 + (z - z^ 2 

where (x,y,z) are the coordinates of the reflection point O r . The 
reflected orthonormal unit base vectors are given by 

/■S 1 ^ A 

x^ = x^ - 2(xj • n) n (C.2a) 

x 2 = x 2 - 2(x^ • n) n (C.2b) 

~ x(x p - x) + y(y~ - y) + z(z p - z) 

x!j = - - -- --- — — -- --- - - (C.2c) 

]/ ((x 2 - x) 2 + (y 2 - y) 2 + (z 2 - z) 2 ) 

The vectors (C.2a) to (c) are chosen to satisfy a left-hand system, 

~r ~ r ~ r 
i .e . , x^ x x^ = -Xy 


This selection is invariant with respect to the 
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.r „r , 


reflected curvatures <R^ , R£>. The three orthonormal base vectors of 
the reflector surface at O r are chosen to be 


x + z f 

x s * 

X 1 " 


J <1 * 


(C.3a) 


x s x x s 
~s _ x 3 X 1 

X'i — 


x 3 x X 1 


(C.3b) 


where 


;« . ; . y * z 


yy 


* fy * » 


f _ d f(x,y) 
x ~ dx 

f d f(x,y) 
y = dy 

z = f(x,y) reflectro surface function 
Let define the following parameters be defined, 


(C.3c) 


(C. 3d) 

(C. 3e) 
(C.3f) 


p’ = X I • x n m - 1.2,3, n=l ,2,3. (C.4) 

mn m n 

Substituting Eqs. < C . la) to <c) and (C.3a) to (c) into (C.4) we obtain. 


(z - z 1 ) - f x ( x - x 1 ) 


1 1 


vAn + fb s/ 


((x - x^ ) 2 + (z - z^ ) 2 ) 


(C.5a) 


-f y (x - x ^ ) + f x (z - z^) 


12 


y 


A(1 + f 2 ) J [(X - X] ) 2 + (Z - 2-j ) 2 ] 


(C. 5b) 


21 


~(x - x 1 ) ( y - - f ^ ( y - y^)(i - 


/ T /f 2 2 2 2*2 2 Z\ 

V ( 1 ♦ M V\(x ~ (y - y 1 ) + C ( Z - z,> + (x - Xj) ] + (y - y ^ ) (z - z^ J 


(C.5c) 



REFLECTOR 



z 


C.l Principal vector directions of the incident and 
reflected rays 
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22 


f f (x - x )(y - y ) +■ (1 + f ) [(z - z ) 2 ♦ (x - x ) 2 ] - f (y - y )(z - z ) 

X V j ] X 1 1 ^ ] ] 

\/5777) \/v x " > 2 (y - y, ) 2 + C(z - ) 2 + (x - x 1 ) 2 ] 2 + (y - y 1 ) 2 (z - z i ) 2 ^ 

(C.5d) 


.1 1 


33 = dl A<f x <x " V + f y (y y l } <z " Z 1 


)) 


(C.5e) 


where 


. ,-2 .2 , N — 1 / 2 

A - (f x + f y + n 


(C.5f) 


Let the 2x2 matrix P be defined for the first four parameters in 
Eq. (C.5a) to (d). If we define another set of parameters by the 
the following definition, 

PL = *1 • m= 1 ,2,3 and n-1,2,3 (C.6) 

mn m n 

it can be shown that 


p r _ p’ 

K 11 " 'll 


p r _ p 1 

12 " 12 


p r _ p 1 

r 21 r 21 


P r - P 1 
r 22 *22 


(C.7a) 

(C.7b) 

(C.7c) 

(C.7d) 


The curvature matrix of incident ray bundle is given by 


Q 1 = 


dj 0 


T'j 


(C.8) 


The curvature matrix of the reflector surface at the reflecting point 
is given by 
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[eA 2 (G - F 2 /E) A(eF - f E ) / E 

Q s _ (C.9) 

U(eF - fE)/E A 2 (gE - 2fF + eF 2 /E)J 

where 


= 1 + f 2 , 

F = f f , G = 

1 + f 2 


x y ’ 

y 

e = - Af 

. f = - Af , 

g = Af 
y yy 

XX 


d 2 f(x,y) 

f d 2 f ( x , v) 

r. d 2 f(x,v) 

‘ d* 2 

r xy " dx dy 

V' dy 2 


The desired reflected curvature matrix Q r may be computed from the 
fol lowing equation 

Q r - Q 1 +2 P^aP 1 ) 1 )" 1 0 S (P 1 )“ 1 (C.10) 

in general Q r is a 2x2 matrix, with elements denoted by Q 1 1 , 

Q 1 2 * ^21 Q 22 • 

Then radii of curvature of the reflected wavefront is given by the 
fol lowing expression, 



(C.ll) 
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APPENDIX D 

FOCAL PLANE FIELDS COMPUTATION BY THE RECIPROCITY PRINCIPLE 
It is required to calculate the focal plane electric field at 
pont 1 produced by an incident uniform plane wave (Fig. (D-l) from an 
observation direction ^gEAM’^BEAM^ ' ** assume that the uniform 

plane wave is generated by a dipole source located at point 2 in the 
far field of the distorted reflector antenna. The dipole source at 
point 2 has the fol lowing current distribution, 

3 2 (r> = u 2 5 (r - r 2 > (D.l) 

where u 2 is a unit vector describing the polarization of the dipole. 
In the neighborhood of any distorted reflector surface point the 
incident field is a local uniform plane wave. This uniform plane wave 
can be described by, 

E 1 = u 2 e +;jK " r (D.2) 

where 

■¥ ~ ~ ~ 
r = xx + yy + zz 

(x,y,z) any point on the distorted reflector surface 

^ A ^ 

K - K (sin Q cos <j> x + sin e s i n 4> y + cos 0 z) 

9 = 0 BEAM and * = ^BEAM 
K = 2 tt/X 

the received focal plane field at point 1 is denoted by. 
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OBSERVATION 

DIRECTION 



(0 BEAM' <l) BEAM ) 


TRANSMITTING 
DIPOLE J2 


D.l Uniform plane wave Incident on the reflector surface 


«<\ 
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(D. 3) 


v p 


(fy. 


*/ P <?, 


Let us now assume that the distorted reflector is being 
considered in the transmitting mode with a dipole source J] located 
at point 1 in the focal plane of the reflector (Fig. D.2). The dipole 
source at point 1 have the following current distribution 

^(r) = u 1 S (r - rj) (D.4) 

where U] is a unit vector describing the transmitting polarization 
of the dipole source In the focal plane. With this dipole as the 
radiating element, the distorted antenna produces a far zone electric 
field at point 2 denoted by 

f 2 tran <? 2 >. a 2 tran <? 2 > (d.5> 


The Lorentz reciprocity theorem relates the transmitting fields 
on Eqs. (D.3) and (D.5) with the respective sources described by Eqs. 
(D.l) and (D.4) respectively. Mathematically the reciprocity principle 
can be expressed as follows, 


'?, fP . a 2 trans - ? 2 trans . a, fp ] . os 


[?, fp • \ - t 2 trans . j, ♦ a 2 trans . a, . a, fp . a 2 ] dv <o.6> 

J J V 

V 

for the particular case under consideration one can use the following 
assumptions. 


^ - 0 and ^ fP 



x r 


Z 


3 tran 

"2 



x r 


(D. 7) 
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U2: (RECEIVING DIPOLE) 



/ 


D.2 Illumination of the reflector surface by a dipole source on 
the focal plane 
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under this assumptions the following result is obtained. 

0 = [[[ [t, fP • H 2 - ? 2 tran . 3,] dv <D.8> 

V 

substituting Eq. (D.l) and (D.4) into Eq. (D.8) yields the final result, 

?, fP • G 2 - ? 2 tran • J, (D.9) 

This result shows that by reciprocity, the focal plane fields 
produced by an Incident uniform plane wave on the distorted reflector 
surface can be obtained indirectly from the far zone electric field 
generated by a dipole source at a specified location in the focal 
plane of the distorted reflector. 



Report Documentation Page 

Space Administration 

1 . Report No. 2. Government Accession No. 

NASA TM- 100826 

3. Recipient’s Catalog No. 

4. Title and Subtitle 

Active Feed Array Compensation for Reflector 
Antenna Surface Distortions 

5. Report Date 

June 1988 

6. Performing Organization Code 

7. Author(s) 

Roberto J. Acosta 

8. Performing Organization Report No. 

E— 401 5 

10. Work Unit No. 

650-60-20 

9. Performing Organization Name and Address 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 

1 1 . Contract or Grant No. 

13. Type of Report and Period Covered 

Technical Memorandum 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546-0001 

14. Sponsoring Agency Code 

15. Supplementary Notes 


This report is a dissertation presented to The University of Akron in partial 
fulfillment of the requirements for the degree Doctor of Philosophy in May 1988. 


16. Abstract 

The feasibility of electromagnetic compensation for reflector antenna surface 
distortions (non-random) is investigated. The performance characteristics (gain, 
sidelobe levels, etc.) of large satellite communication reflector antenna systems 
degrade as the reflector surface distorts, mainly due to thermal effects from a 
solar radiation. The technique developed in this dissertation can be used to 
maintain the antenna boresight directivity and sidelobe level independent of 
thermal effects on the reflector surface. With the advent of monolithic micro- 
wave integrated circuits (MMIC), a greater flexibility in array-fed reflector 
antenna systems can be achieved. MMIC arrays provide independent control of 
amplitude and phase for each of the many radiating elements in the feed array. 

By assuming a known surface distortion profile, a simulation study is carried 
out to examine the antenna performance as a function of feed array size and 
number of elements. Results indicate that the compensation technique can effec- 
tively control boresight directivity (within 1-3 dB of the undistorted) and side- 
lobe level (within 1-5 dB of the undistorted) under peak surface distortion in 
the order of tenth of a wavelength. Furthermore, the computational ease of the 
compensation technique allows it to be implemented adaptively in a large satel- 


lite reflector antenna system. 

1 7. Key Words (Suggested by Author(s)) 

Compensation; Distorted antenna; 
Radiation pattern; Feed array; Conjugate 
field matching 

18. Distribution Statement 

Unclassified - Unlimited 
Subject Category 32 

19. Security Classif. (of this report) 

Unclassified 

20. Security Classif. (of this page) 

Unclassified 

21. No of pages 

178 

22. Price* 

A0 9 


NASA FORM 1626 OCT 86 


*For sale by the National Technical Information Service, Springfield, Virginia 22161 




















National Aeronautics and second class mail 

Space Administration 

Lewis Research Center address correction requested 

Cleveland, Ohio 44135 

Official Business 
Penalty far Private Use $300 



Postage and Fees Paid 
National Aeronautics and 
Space Administration 
NASA-451 





